Rhizotron
dframe_rhiz = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/Rizhotron.txt")
dframe_rhiz_ufscar = dframe_rhiz %>%
filter(Intituto=='UFSCar') %>% select(-Manejo,-Tratamentos) %>%
mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')),
Rhizotron_2022=10*Rhizotron_2022) %>%
pivot_longer(cols = c("Rhizotron_2021","Rhizotron_2022"),names_to = 'Safra',
values_to = 'Rhizotron') %>%
mutate(Safra=gsub('Rhizotron_','',Safra))
dframe_rhiz_ccsm = dframe_rhiz %>%
filter(Intituto=='CCSM') %>% select(-Tratamento,-Tratamentos) %>%
mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV'))) %>%
pivot_longer(cols = c("Rhizotron_2021","Rhizotron_2022"),names_to = 'Safra',
values_to = 'Rhizotron') %>%
mutate(Safra=gsub('Rhizotron_','',Safra))
UFSCar
dframe_summ_ufscar = dframe_rhiz_ufscar %>%
select(Safra,RegHidrico,Tratamento,Rhizotron) %>%
group_by(Safra,RegHidrico,Tratamento) %>%
summarise(Mean = mean(Rhizotron,na.rm=T),
SD = sd(Rhizotron,na.rm=T)) %>%
mutate(`Mean ± SD`=paste(round(Mean,5),'±',round(SD,5)))
Graficos
Gráfico De Barras
Option 1
yname = 'Média de rhizotron'
dframe_summ_ufscar %>% ggplot(
aes(x=Tratamento,
fill=RegHidrico,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
facet_wrap(~Safra)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)

Option 2
dframe_summ_ufscar %>% ggplot(
aes(x=RegHidrico,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
facet_wrap(~Safra)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)

Gráfico de Pontos + Média e Erro
Option 1
dframe_rhiz_ufscar %>% ggplot()+
facet_wrap(~Safra,scales='free')+
geom_point(aes(x=Tratamento,y=Rhizotron,
color=RegHidrico,
shape=RegHidrico),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(x=Tratamento,y=Mean,
color=RegHidrico,
shape=RegHidrico),
data=dframe_summ_ufscar,
position = position_dodge(width = .75))+
geom_errorbar(aes(x=Tratamento,y=Mean,
color=RegHidrico,
linetype=RegHidrico,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ_ufscar,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()

Option 2
dframe_rhiz_ufscar %>% ggplot()+
facet_wrap(~Safra,scales='free')+
geom_point(aes(x=RegHidrico,y=Rhizotron,
color=Tratamento,
shape=Tratamento),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(x=RegHidrico,y=Mean,
color=Tratamento,
shape=Tratamento),
data=dframe_summ_ufscar,
position = position_dodge(width = .75))+
geom_errorbar(aes(x=RegHidrico,y=Mean,
color=Tratamento,
linetype=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ_ufscar,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()

Tabela Resumo
dframe_summ_ufscar %>% select(Safra,RegHidrico,Tratamento,`Mean ± SD`) %>% kable_cust(cap = 'Média e Desvio Padrao - Rhizotron')
Média e Desvio Padrao - Rhizotron
|
Safra
|
RegHidrico
|
Tratamento
|
Mean ± SD
|
|
2021
|
S
|
BRF_PF+FD
|
0.01001 ± 0.00631
|
|
2021
|
S
|
H10+FD
|
0.0105 ± 0.00623
|
|
2021
|
S
|
H10+H152
|
0.03834 ± 0.01675
|
|
2021
|
TRF
|
BRF_PF+FD
|
0.02045 ± 0.02714
|
|
2021
|
TRF
|
H10+FD
|
0.00918 ± 0.01088
|
|
2021
|
TRF
|
H10+H152
|
0.01859 ± 0.0209
|
|
2021
|
TRV
|
BRF_PF+FD
|
0.00808 ± 0.00513
|
|
2021
|
TRV
|
H10+FD
|
0.01402 ± 0.0136
|
|
2021
|
TRV
|
H10+H152
|
0.02275 ± 0.00243
|
|
2022
|
S
|
BRF_PF+FD
|
0.03457 ± 0.00308
|
|
2022
|
S
|
H10+FD
|
0.03134 ± 0.00449
|
|
2022
|
S
|
H10+H152
|
0.04866 ± 0.00287
|
|
2022
|
TRF
|
BRF_PF+FD
|
0.02024 ± 0.00499
|
|
2022
|
TRF
|
H10+FD
|
0.01622 ± 0.00192
|
|
2022
|
TRF
|
H10+H152
|
0.02094 ± 0.00723
|
|
2022
|
TRV
|
BRF_PF+FD
|
0.01795 ± 0.00295
|
|
2022
|
TRV
|
H10+FD
|
0.02425 ± 0.00266
|
|
2022
|
TRV
|
H10+H152
|
0.02961 ± 0.00519
|
Modelo
2021
df_tk <- dframe_rhiz_ufscar %>% filter(Safra=='2021')
fit_lm <- lm(data=df_tk,
formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)
fit_aov <- aov(data=df_tk,
formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)
summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)
summ_tk <- df_tk %>% group_by(RegHidrico,Tratamento) %>%
summarise(mean=mean(Rhizotron)) %>% arrange(-mean)
summ_tk$Group =as.data.frame.list(mtcomp$`RegHidrico:Tratamento`)[,1]
summ_tk %>% kable_cust()
|
RegHidrico
|
Tratamento
|
mean
|
Group
|
|
S
|
H10+H152
|
0.03834
|
a
|
|
TRV
|
H10+H152
|
0.02275
|
a
|
|
TRF
|
BRF_PF+FD
|
0.02045
|
a
|
|
TRF
|
H10+H152
|
0.01859
|
a
|
|
TRV
|
H10+FD
|
0.01402
|
a
|
|
S
|
H10+FD
|
0.01050
|
a
|
|
S
|
BRF_PF+FD
|
0.01001
|
a
|
|
TRF
|
H10+FD
|
0.00918
|
a
|
|
TRV
|
BRF_PF+FD
|
0.00808
|
a
|
#tx <- with(data = df_tk,interaction(RegHidrico,Tratamento))
#amod <- aov(Rhizotron ~ RegHidrico+Tratamento+tx, data=df_tk)
#Tukey = agricolae::HSD.test(amod,trt = 'tx',group=TRUE)
#Tukey$groups %>% kable_cust()
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
summ_aov
## Df Sum Sq Mean Sq F value Pr(>F)
## RegHidrico 2 0.000107 0.0000534 0.257 0.7761
## Tratamento 2 0.001277 0.0006387 3.077 0.0709 .
## RegHidrico:Tratamento 4 0.000847 0.0002116 1.020 0.4236
## Residuals 18 0.003736 0.0002076
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
0.01001
|
0.00832
|
1.20307
|
0.24454
|
|
RegHidricoTRF
|
0.01044
|
0.01176
|
0.88738
|
0.38658
|
|
RegHidricoTRV
|
-0.00192
|
0.01176
|
-0.16355
|
0.87191
|
|
TratamentoH10+FD
|
0.00049
|
0.01176
|
0.04153
|
0.96733
|
|
TratamentoH10+H152
|
0.02833
|
0.01176
|
2.40859
|
0.02695
|
|
RegHidricoTRF:TratamentoH10+FD
|
-0.01176
|
0.01664
|
-0.70673
|
0.48878
|
|
RegHidricoTRV:TratamentoH10+FD
|
0.00544
|
0.01664
|
0.32723
|
0.74727
|
|
RegHidricoTRF:TratamentoH10+H152
|
-0.03019
|
0.01664
|
-1.81452
|
0.08630
|
|
RegHidricoTRV:TratamentoH10+H152
|
-0.01367
|
0.01664
|
-0.82149
|
0.42211
|
2022
df_tk <- dframe_rhiz_ufscar %>% filter(Safra=='2022')
fit_lm <- lm(data=df_tk,
formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)
fit_aov <- aov(data=df_tk,
formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)
summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)
summ_tk <- df_tk %>% group_by(RegHidrico,Tratamento) %>%
summarise(mean=mean(Rhizotron)) %>% arrange(-mean)
summ_tk$Group = as.data.frame.list(mtcomp$`RegHidrico:Tratamento`)[,1]
summ_tk %>% kable_cust()
|
RegHidrico
|
Tratamento
|
mean
|
Group
|
|
S
|
H10+H152
|
0.04866
|
a
|
|
S
|
BRF_PF+FD
|
0.03457
|
b
|
|
S
|
H10+FD
|
0.03134
|
bc
|
|
TRV
|
H10+H152
|
0.02961
|
bcd
|
|
TRV
|
H10+FD
|
0.02425
|
bcde
|
|
TRF
|
H10+H152
|
0.02094
|
cde
|
|
TRF
|
BRF_PF+FD
|
0.02024
|
cde
|
|
TRV
|
BRF_PF+FD
|
0.01795
|
de
|
|
TRF
|
H10+FD
|
0.01622
|
e
|
#tx <- with(data = df_tk,interaction(RegHidrico,Tratamento))
#amod <- aov(Rhizotron ~ RegHidrico+Tratamento+tx, data=df_tk)
#
#Tukey = agricolae::HSD.test(amod,trt = 'tx',group=TRUE)
#Tukey$groups %>% kable_cust()
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
summ_aov
## Df Sum Sq Mean Sq F value Pr(>F)
## RegHidrico 2 0.0017678 0.0008839 49.258 5.01e-08 ***
## Tratamento 2 0.0004839 0.0002419 13.483 0.000264 ***
## RegHidrico:Tratamento 4 0.0002684 0.0000671 3.739 0.021999 *
## Residuals 18 0.0003230 0.0000179
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
0.03457
|
0.00245
|
14.13354
|
0.00000
|
|
RegHidricoTRF
|
-0.01433
|
0.00346
|
-4.14327
|
0.00061
|
|
RegHidricoTRV
|
-0.01661
|
0.00346
|
-4.80345
|
0.00014
|
|
TratamentoH10+FD
|
-0.00323
|
0.00346
|
-0.93337
|
0.36298
|
|
TratamentoH10+H152
|
0.01409
|
0.00346
|
4.07497
|
0.00071
|
|
RegHidricoTRF:TratamentoH10+FD
|
-0.00079
|
0.00489
|
-0.16097
|
0.87391
|
|
RegHidricoTRV:TratamentoH10+FD
|
0.00953
|
0.00489
|
1.94779
|
0.06721
|
|
RegHidricoTRF:TratamentoH10+H152
|
-0.01339
|
0.00489
|
-2.73656
|
0.01356
|
|
RegHidricoTRV:TratamentoH10+H152
|
-0.00244
|
0.00489
|
-0.49902
|
0.62381
|
CCSM
dframe_summ_ccsm = dframe_rhiz_ccsm %>%
select(Safra,Manejo,Copa,PortaEnxerto,Rhizotron) %>%
group_by(Safra,Manejo,Copa,PortaEnxerto) %>%
summarise(Mean = mean(Rhizotron,na.rm=T),
SD = sd(Rhizotron,na.rm=T)) %>%
mutate(`Mean ± SD`=paste(round(Mean,5),'±',round(SD,5)))
Graficos
Gráfico de Barras
Option 1
yname = 'Média de rhizotron'
dframe_summ_ccsm %>% ggplot(
aes(x=PortaEnxerto,
fill=Copa,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
facet_wrap(Safra~Manejo,nrow=1)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)

Option 2
dframe_summ_ccsm %>% ggplot(
aes(x=Copa,
fill=PortaEnxerto,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
facet_wrap(Safra~Manejo,nrow=1)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)

Gráfico de Pontos + Média e Erro
Option 1
dframe_rhiz_ccsm %>% ggplot()+
facet_wrap(~Safra)+
geom_point(aes(x=PortaEnxerto,
y=Rhizotron,
color=Copa,
shape=Manejo),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(x=PortaEnxerto,
y=Mean,
color=Copa,
shape=Manejo),
data=dframe_summ_ccsm,
position = position_dodge(width = .75))+
geom_errorbar(aes(x=PortaEnxerto,
y=Mean,
color=Copa,
shape=Manejo,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ_ccsm,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()

Option 2
dframe_rhiz_ccsm %>% ggplot()+
facet_wrap(~Safra)+
geom_point(aes(x=Copa,
y=Rhizotron,
color=PortaEnxerto,
shape=Manejo),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(x=Copa,
y=Mean,
color=PortaEnxerto,
shape=Manejo),
data=dframe_summ_ccsm,
position = position_dodge(width = .75))+
geom_errorbar(aes(x=Copa,
y=Mean,
color=PortaEnxerto,
shape=Manejo,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ_ccsm,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()

Tabela Resumo
dframe_summ_ccsm %>% select(Safra,Manejo,PortaEnxerto,Copa,`Mean ± SD`) %>% kable_cust(cap = 'Média e Desvio Padrao - Rhizotron')
Média e Desvio Padrao - Rhizotron
|
Safra
|
Manejo
|
PortaEnxerto
|
Copa
|
Mean ± SD
|
|
2021
|
Convencional
|
FD
|
BRF_PF
|
0.00521 ± 0.00401
|
|
2021
|
Convencional
|
H152
|
BRF_PF
|
0.00837 ± 0.0044
|
|
2021
|
Convencional
|
FD
|
IAC10
|
0.0035 ± 0.00158
|
|
2021
|
Convencional
|
H152
|
IAC10
|
0.01423 ± 0.00303
|
|
2021
|
Convencional
|
FD
|
IAC5
|
0.01404 ± 0.01049
|
|
2021
|
Convencional
|
H152
|
IAC5
|
0.00987 ± 0.00401
|
|
2021
|
Ecologico
|
FD
|
BRF_PF
|
0.0088 ± 0.0014
|
|
2021
|
Ecologico
|
H152
|
BRF_PF
|
0.01557 ± 0.00379
|
|
2021
|
Ecologico
|
FD
|
IAC10
|
0.00549 ± 0.00033
|
|
2021
|
Ecologico
|
H152
|
IAC10
|
0.03573 ± 0.01341
|
|
2021
|
Ecologico
|
FD
|
IAC5
|
0.01022 ± 0.01135
|
|
2021
|
Ecologico
|
H152
|
IAC5
|
0.01021 ± 0.00289
|
|
2022
|
Convencional
|
FD
|
BRF_PF
|
0.01011 ± 0.00569
|
|
2022
|
Convencional
|
H152
|
BRF_PF
|
0.02145 ± 0.00038
|
|
2022
|
Convencional
|
FD
|
IAC10
|
0.00386 ± 0.00174
|
|
2022
|
Convencional
|
H152
|
IAC10
|
0.00909 ± 0.00083
|
|
2022
|
Convencional
|
FD
|
IAC5
|
0.00171 ± 6e-05
|
|
2022
|
Convencional
|
H152
|
IAC5
|
0.01492 ± 0.00074
|
|
2022
|
Ecologico
|
FD
|
BRF_PF
|
0.00425 ± 0.00071
|
|
2022
|
Ecologico
|
H152
|
BRF_PF
|
0.02376 ± 0.00298
|
|
2022
|
Ecologico
|
FD
|
IAC10
|
0.00566 ± 0.00022
|
|
2022
|
Ecologico
|
H152
|
IAC10
|
0.00891 ± 0.00142
|
|
2022
|
Ecologico
|
FD
|
IAC5
|
0.00459 ± 0.00013
|
|
2022
|
Ecologico
|
H152
|
IAC5
|
0.01106 ± 0.00222
|
Modelo
2021
df_tk <- dframe_rhiz_ccsm %>% filter(Safra=='2021')
fit_lm <- lm(data=df_tk,
formula = Rhizotron~Manejo+Copa +PortaEnxerto+
Manejo*Copa+Manejo*PortaEnxerto +
Copa*PortaEnxerto+
Manejo*Copa*PortaEnxerto
)
fit_aov <- aov(data=df_tk,
formula = Rhizotron~Manejo+Copa +PortaEnxerto+
Manejo*Copa+#Manejo*PortaEnxerto +
Copa*PortaEnxerto+#+
Manejo*Copa*PortaEnxerto
)
summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)
summ_tk <- df_tk %>% group_by(Manejo,Copa,PortaEnxerto) %>%
summarise(mean=mean(Rhizotron)) %>% arrange(-mean)
summ_tk$Group =as.data.frame.list(mtcomp$`Manejo:Copa:PortaEnxerto`)[,1]
summ_tk %>% kable_cust()
|
Manejo
|
Copa
|
PortaEnxerto
|
mean
|
Group
|
|
Ecologico
|
IAC10
|
H152
|
0.03573
|
a
|
|
Ecologico
|
BRF_PF
|
H152
|
0.01557
|
b
|
|
Convencional
|
IAC10
|
H152
|
0.01423
|
b
|
|
Convencional
|
IAC5
|
FD
|
0.01404
|
b
|
|
Ecologico
|
IAC5
|
FD
|
0.01022
|
b
|
|
Ecologico
|
IAC5
|
H152
|
0.01021
|
b
|
|
Convencional
|
IAC5
|
H152
|
0.00987
|
b
|
|
Ecologico
|
BRF_PF
|
FD
|
0.00880
|
b
|
|
Convencional
|
BRF_PF
|
H152
|
0.00837
|
b
|
|
Ecologico
|
IAC10
|
FD
|
0.00549
|
b
|
|
Convencional
|
BRF_PF
|
FD
|
0.00521
|
b
|
|
Convencional
|
IAC10
|
FD
|
0.00350
|
b
|
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
summ_aov
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0002372 0.0002372 5.620 0.026137 *
## Copa 2 0.0001736 0.0000868 2.056 0.149916
## PortaEnxerto 1 0.0005457 0.0005457 12.928 0.001454 **
## Manejo:Copa 2 0.0002731 0.0001366 3.236 0.057000 .
## Copa:PortaEnxerto 2 0.0007999 0.0004000 9.476 0.000926 ***
## Manejo:PortaEnxerto 1 0.0001861 0.0001861 4.409 0.046447 *
## Manejo:Copa:PortaEnxerto 2 0.0001221 0.0000611 1.447 0.255160
## Residuals 24 0.0010130 0.0000422
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
0.00521
|
0.00375
|
1.38946
|
0.17745
|
|
ManejoEcologico
|
0.00359
|
0.00530
|
0.67625
|
0.50535
|
|
CopaIAC10
|
-0.00171
|
0.00530
|
-0.32305
|
0.74946
|
|
CopaIAC5
|
0.00883
|
0.00530
|
1.66454
|
0.10901
|
|
PortaEnxertoH152
|
0.00316
|
0.00530
|
0.59530
|
0.55722
|
|
ManejoEcologico:CopaIAC10
|
-0.00160
|
0.00750
|
-0.21266
|
0.83339
|
|
ManejoEcologico:CopaIAC5
|
-0.00741
|
0.00750
|
-0.98735
|
0.33332
|
|
ManejoEcologico:PortaEnxertoH152
|
0.00362
|
0.00750
|
0.48201
|
0.63416
|
|
CopaIAC10:PortaEnxertoH152
|
0.00757
|
0.00750
|
1.00905
|
0.32302
|
|
CopaIAC5:PortaEnxertoH152
|
-0.00733
|
0.00750
|
-0.97650
|
0.33856
|
|
ManejoEcologico:CopaIAC10:PortaEnxertoH152
|
0.01589
|
0.01061
|
1.49799
|
0.14717
|
|
ManejoEcologico:CopaIAC5:PortaEnxertoH152
|
0.00054
|
0.01061
|
0.05116
|
0.95962
|
2022
df_tk <- dframe_rhiz_ccsm %>% filter(Safra=='2022')
fit_lm <- lm(data=df_tk,
formula = Rhizotron~
#Manejo+
Copa +PortaEnxerto+
#Manejo*Copa+Manejo*PortaEnxerto +
Copa*PortaEnxerto#+
#Manejo*Copa*PortaEnxerto
)
fit_aov <- aov(data=df_tk,
formula = Rhizotron~
Manejo+Copa +PortaEnxerto+
Manejo*Copa+Manejo*PortaEnxerto +
Copa*PortaEnxerto+
Manejo*Copa*PortaEnxerto
)
summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)
summ_tk <- df_tk %>% group_by(Copa,PortaEnxerto) %>%
summarise(mean=mean(Rhizotron)) %>% arrange(-mean)
summ_tk$Group =as.data.frame.list(mtcomp$`Copa:PortaEnxerto`)[,1]
summ_tk %>% kable_cust()
|
Copa
|
PortaEnxerto
|
mean
|
Group
|
|
BRF_PF
|
H152
|
0.02261
|
a
|
|
IAC5
|
H152
|
0.01299
|
b
|
|
IAC10
|
H152
|
0.00900
|
c
|
|
BRF_PF
|
FD
|
0.00718
|
cd
|
|
IAC10
|
FD
|
0.00476
|
de
|
|
IAC5
|
FD
|
0.00315
|
e
|
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
par(mfrow=c(2,2))
plot(fit_aov)

par(mfrow=c(1,1))
summ_aov
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0000021 0.0000021 0.480 0.495002
## Copa 2 0.0004488 0.0002244 50.641 2.44e-09 ***
## PortaEnxerto 1 0.0008701 0.0008701 196.364 4.73e-13 ***
## Manejo:Copa 2 0.0000099 0.0000050 1.122 0.342089
## Manejo:PortaEnxerto 1 0.0000001 0.0000001 0.017 0.896199
## Copa:PortaEnxerto 2 0.0001877 0.0000938 21.174 5.02e-06 ***
## Manejo:Copa:PortaEnxerto 2 0.0000871 0.0000435 9.827 0.000763 ***
## Residuals 24 0.0001063 0.0000044
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
0.00718
|
0.00107
|
6.71937
|
0.00000
|
|
CopaIAC10
|
-0.00242
|
0.00151
|
-1.60200
|
0.11963
|
|
CopaIAC5
|
-0.00403
|
0.00151
|
-2.66480
|
0.01228
|
|
PortaEnxertoH152
|
0.01543
|
0.00151
|
10.20594
|
0.00000
|
|
CopaIAC10:PortaEnxertoH152
|
-0.01118
|
0.00214
|
-5.23293
|
0.00001
|
|
CopaIAC5:PortaEnxertoH152
|
-0.00559
|
0.00214
|
-2.61646
|
0.01378
|
Qualidade do fruto
dframe_qf = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/QualidadeFruto.txt")
dframe_qf_ufscar = dframe_qf %>%
filter(Instituto=='UFSCar') %>% select(-Manejo) %>%
mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))
dframe_qf_ccsm = dframe_qf %>%
filter(Instituto=='CCSM') %>% select(-Tratamento) %>%
mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))
summ_qf <- function(dframe,variable_name = ''){
if (variable_name==''){
stop()
}
dframe = dframe %>% select(one_of(variable_name),RegHidrico,Tratamento)
colnames(dframe)[1]<-'Y'
out<- list()
dframe_summ <- dframe %>% group_by(RegHidrico,Tratamento) %>%
summarise(Mean = mean(Y,na.rm=T),
SD = sd(Y,na.rm=T)) %>%
mutate(`Mean ± SD`=paste(round(Mean,2),'±',round(SD,2)))
out$dframe_summ <- dframe_summ
out$barp <- dframe_summ %>% ggplot(
aes(x=RegHidrico,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_2 <- dframe_summ %>% ggplot(
aes(x=Tratamento,
fill=RegHidrico,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_error <- dframe_summ %>% ggplot(
aes(x=RegHidrico,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
geom_errorbar(
aes(x=RegHidrico,y=Mean,
color=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_error_2 <- dframe_summ %>% ggplot(
aes(x=Tratamento,
fill=RegHidrico,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
geom_errorbar(
aes(x=Tratamento,y=Mean,
color=RegHidrico,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$mean_desv_plot <-
dframe %>% ggplot()+
geom_point(aes(
x=RegHidrico,
y=Y,
color=Tratamento,
shape=Tratamento),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(
x=RegHidrico,
y=Mean,
color=Tratamento,
shape=Tratamento),
data=dframe_summ,
position = position_dodge(width = .75))+
geom_errorbar(aes(
x=RegHidrico,
y=Mean,
color=Tratamento,
shape=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$mean_desv_plot_2 <-
dframe %>% ggplot()+
geom_point(aes(
x=Tratamento,
y=Y,
color=RegHidrico,
shape=RegHidrico),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(
x=Tratamento,
y=Mean,
color=RegHidrico,
shape=RegHidrico),
data=dframe_summ,
position = position_dodge(width = .75))+
geom_errorbar(aes(
x=Tratamento,
y=Mean,
color=RegHidrico,
shape=RegHidrico,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
return(out)
}
summ_qf_ccsm <- function(dframe,variable_name = ''){
if (variable_name==''){
stop()
}
dframe = dframe %>% select(one_of(variable_name),Manejo,Copa,PortaEnxerto)
colnames(dframe)[1]<-'Y'
out<- list()
dframe = dframe %>% mutate(Tratamento = paste(Copa,'+',PortaEnxerto,sep=''))
dframe_summ <- dframe %>% group_by(Manejo,Tratamento) %>%
summarise(Mean = mean(Y,na.rm=T),
SD = sd(Y,na.rm=T)) %>%
mutate(`Mean ± SD`=paste(round(Mean,2),'±',round(SD,2)))
out$dframe_summ <- dframe_summ
out$barp <- dframe_summ %>% ggplot(
aes(x=Manejo,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_2 <- dframe_summ %>% ggplot(
aes(x=Tratamento,
fill=Manejo,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_error <- dframe_summ %>% ggplot(
aes(x=Manejo,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
geom_errorbar(
aes(x=Manejo,y=Mean,
color=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_error_2 <- dframe_summ %>% ggplot(
aes(x=Tratamento,
fill=Manejo,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
geom_errorbar(
aes(x=Tratamento,y=Mean,
color=Manejo,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$mean_desv_plot <-
dframe %>% ggplot()+
geom_point(aes(
x=Manejo,
y=Y,
color=Tratamento,
shape=Tratamento),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(
x=Manejo,
y=Mean,
color=Tratamento,
shape=Tratamento),
data=dframe_summ,
position = position_dodge(width = .75))+
geom_errorbar(aes(
x=Manejo,
y=Mean,
color=Tratamento,
shape=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$mean_desv_plot_2 <-
dframe %>% ggplot()+
geom_point(aes(
x=Tratamento,
y=Y,
color=Manejo,
shape=Manejo),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(
x=Tratamento,
y=Mean,
color=Manejo,
shape=Manejo),
data=dframe_summ,
position = position_dodge(width = .75))+
geom_errorbar(aes(
x=Tratamento,
y=Mean,
color=Manejo,
shape=Manejo,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
return(out)
}
### model funcitons =
dframe_aux = dframe_qf_ccsm
response='Massa'
fit_qf_ccsm <- function(dframe_aux,response){
dframe_aux = dframe_aux %>% select(Manejo,PortaEnxerto,Copa,one_of(response))
out_models = list()
out_models$models <- list()
## Full Model
out_models$models$fit_full = lm(data=dframe_aux,get(response)~Manejo*PortaEnxerto*Copa)
out_models$models$fit_full_aov = aov(data=dframe_aux,get(response)~Manejo*PortaEnxerto*Copa)
## All 2*2 ints
out_models$models$fit_noInt3 = lm(data=dframe_aux,get(response) ~ Manejo*Copa +Manejo*PortaEnxerto + PortaEnxerto*Copa)
out_models$models$fit_noInt3_aov = aov(data=dframe_aux,get(response)~ Manejo*Copa +Manejo*PortaEnxerto + PortaEnxerto*Copa)
## Only PortaEnxerto*Copa
out_models$models$fit_CopPort = lm(data=dframe_aux,get(response)~ Manejo + PortaEnxerto*Copa)
out_models$models$fit_CopPort_aov = aov(data=dframe_aux,get(response)~ Manejo + PortaEnxerto*Copa)
## Just main effects
out_models$models$fit_mainEffect = lm(data=dframe_aux,get(response)~ Manejo + PortaEnxerto + Copa)
out_models$models$fit_mainEffect_aov = aov(data=dframe_aux,get(response)~ Manejo + PortaEnxerto + Copa)
### Tukey
out_models$tukey <- list()
tk <- TukeyHSD(out_models$models$fit_full_aov)
mtcomp <- multcompView::multcompLetters4(out_models$models$fit_full_aov,tk)
out_models$tukey$fit_full_aov <-
dframe_aux %>% group_by(Manejo,PortaEnxerto,Copa) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$`Manejo:PortaEnxerto:Copa`$Letters)
tk <- TukeyHSD(out_models$models$fit_noInt3_aov)
mtcomp <- multcompView::multcompLetters4(out_models$models$fit_noInt3_aov,tk)
out_models$tukey$fit_noInt3_aov <- list(
dframe_aux %>% group_by(Manejo,Copa) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$`Manejo:Copa`$Letters),
dframe_aux %>% group_by(Manejo,PortaEnxerto) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$`Manejo:PortaEnxerto`$Letters),
dframe_aux %>% group_by(Copa,PortaEnxerto) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$`Copa:PortaEnxerto`$Letters)
)
tk <- TukeyHSD(out_models$models$fit_CopPort_aov)
mtcomp <- multcompView::multcompLetters4(out_models$models$fit_CopPort_aov,tk)
out_models$tukey$fit_CopPort_aov <- list(
dframe_aux %>% group_by(Manejo,Copa) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$`PortaEnxerto:Copa`$Letters),
dframe_aux %>% group_by(Manejo) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$Manejo$Letters),
dframe_aux %>% group_by(Copa) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$Copa$Letters),
dframe_aux %>% group_by(PortaEnxerto) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$PortaEnxerto$Letters)
)
tk <- TukeyHSD(out_models$models$fit_mainEffect_aov)
mtcomp <- multcompView::multcompLetters4(out_models$models$fit_mainEffect_aov,tk)
out_models$tukey$fit_mainEffect_aov <- list(
dframe_aux %>% group_by(Manejo) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$Manejo$Letters),
dframe_aux %>% group_by(Copa) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$Copa$Letters),
dframe_aux %>% group_by(PortaEnxerto) %>%
summarise(Mean=mean(get(response))) %>%
arrange(-Mean) %>% bind_cols(Group = mtcomp$PortaEnxerto$Letters)
)
return(out_models)
}
UFScar
Massa
results <- summ_qf(dframe_qf_ufscar,variable_name = "Massa")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
81.66667
|
7.09460
|
81.67 ± 7.09
|
|
S
|
FD+H10
|
47.00000
|
2.64575
|
47 ± 2.65
|
|
S
|
H14+H10
|
69.33333
|
3.51188
|
69.33 ± 3.51
|
|
S
|
H152+H10
|
69.33333
|
6.11010
|
69.33 ± 6.11
|
|
S
|
H228+H10
|
47.66667
|
0.57735
|
47.67 ± 0.58
|
|
TRF
|
FD+BRS_PF
|
86.66667
|
8.14453
|
86.67 ± 8.14
|
|
TRF
|
FD+H10
|
82.06667
|
4.10041
|
82.07 ± 4.1
|
|
TRF
|
H14+H10
|
76.00000
|
2.64575
|
76 ± 2.65
|
|
TRF
|
H152+H10
|
79.50000
|
8.50000
|
79.5 ± 8.5
|
|
TRF
|
H228+H10
|
83.00000
|
8.18535
|
83 ± 8.19
|
|
TRV
|
FD+BRS_PF
|
71.66667
|
1.15470
|
71.67 ± 1.15
|
|
TRV
|
FD+H10
|
75.33333
|
1.15470
|
75.33 ± 1.15
|
|
TRV
|
H14+H10
|
54.66667
|
8.50490
|
54.67 ± 8.5
|
|
TRV
|
H152+H10
|
74.66667
|
3.51188
|
74.67 ± 3.51
|
|
TRV
|
H228+H10
|
68.33333
|
7.02377
|
68.33 ± 7.02
|
Altura
results <- summ_qf(dframe_qf_ufscar,variable_name = "Altura")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
5.80000
|
0.20000
|
5.8 ± 0.2
|
|
S
|
FD+H10
|
4.86667
|
0.11547
|
4.87 ± 0.12
|
|
S
|
H14+H10
|
5.53333
|
0.05774
|
5.53 ± 0.06
|
|
S
|
H152+H10
|
5.40000
|
0.20000
|
5.4 ± 0.2
|
|
S
|
H228+H10
|
5.06667
|
0.05774
|
5.07 ± 0.06
|
|
TRF
|
FD+BRS_PF
|
5.86667
|
0.28868
|
5.87 ± 0.29
|
|
TRF
|
FD+H10
|
5.83333
|
0.15275
|
5.83 ± 0.15
|
|
TRF
|
H14+H10
|
5.60000
|
0.00000
|
5.6 ± 0
|
|
TRF
|
H152+H10
|
5.70000
|
0.10000
|
5.7 ± 0.1
|
|
TRF
|
H228+H10
|
5.73333
|
0.25166
|
5.73 ± 0.25
|
|
TRV
|
FD+BRS_PF
|
5.43333
|
0.11547
|
5.43 ± 0.12
|
|
TRV
|
FD+H10
|
5.50000
|
0.10000
|
5.5 ± 0.1
|
|
TRV
|
H14+H10
|
5.13333
|
0.11547
|
5.13 ± 0.12
|
|
TRV
|
H152+H10
|
5.46667
|
0.32146
|
5.47 ± 0.32
|
|
TRV
|
H228+H10
|
5.53333
|
0.28868
|
5.53 ± 0.29
|
Largura
results <- summ_qf(dframe_qf_ufscar,variable_name = "Largura")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
5.26667
|
0.11547
|
5.27 ± 0.12
|
|
S
|
FD+H10
|
4.36667
|
0.05774
|
4.37 ± 0.06
|
|
S
|
H14+H10
|
5.00000
|
0.10000
|
5 ± 0.1
|
|
S
|
H152+H10
|
5.03333
|
0.15275
|
5.03 ± 0.15
|
|
S
|
H228+H10
|
4.40000
|
0.10000
|
4.4 ± 0.1
|
|
TRF
|
FD+BRS_PF
|
5.43333
|
0.15275
|
5.43 ± 0.15
|
|
TRF
|
FD+H10
|
5.33333
|
0.11547
|
5.33 ± 0.12
|
|
TRF
|
H14+H10
|
5.23333
|
0.05774
|
5.23 ± 0.06
|
|
TRF
|
H152+H10
|
5.40000
|
0.10000
|
5.4 ± 0.1
|
|
TRF
|
H228+H10
|
5.30000
|
0.17321
|
5.3 ± 0.17
|
|
TRV
|
FD+BRS_PF
|
5.10000
|
0.00000
|
5.1 ± 0
|
|
TRV
|
FD+H10
|
5.16667
|
0.05774
|
5.17 ± 0.06
|
|
TRV
|
H14+H10
|
4.63333
|
0.30551
|
4.63 ± 0.31
|
|
TRV
|
H152+H10
|
4.93333
|
0.28868
|
4.93 ± 0.29
|
|
TRV
|
H228+H10
|
4.96667
|
0.15275
|
4.97 ± 0.15
|
AL
results <- summ_qf(dframe_qf_ufscar,variable_name = "AL")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
1.10333
|
0.02082
|
1.1 ± 0.02
|
|
S
|
FD+H10
|
1.11667
|
0.02517
|
1.12 ± 0.03
|
|
S
|
H14+H10
|
1.10667
|
0.01155
|
1.11 ± 0.01
|
|
S
|
H152+H10
|
1.07333
|
0.03055
|
1.07 ± 0.03
|
|
S
|
H228+H10
|
1.15333
|
0.04041
|
1.15 ± 0.04
|
|
TRF
|
FD+BRS_PF
|
1.08333
|
0.02517
|
1.08 ± 0.03
|
|
TRF
|
FD+H10
|
1.09333
|
0.02082
|
1.09 ± 0.02
|
|
TRF
|
H14+H10
|
1.07333
|
0.01155
|
1.07 ± 0.01
|
|
TRF
|
H152+H10
|
1.05667
|
0.00577
|
1.06 ± 0.01
|
|
TRF
|
H228+H10
|
1.08333
|
0.02082
|
1.08 ± 0.02
|
|
TRV
|
FD+BRS_PF
|
1.06667
|
0.02309
|
1.07 ± 0.02
|
|
TRV
|
FD+H10
|
1.06667
|
0.02309
|
1.07 ± 0.02
|
|
TRV
|
H14+H10
|
1.11000
|
0.05000
|
1.11 ± 0.05
|
|
TRV
|
H152+H10
|
1.11000
|
0.01000
|
1.11 ± 0.01
|
|
TRV
|
H228+H10
|
1.11333
|
0.03055
|
1.11 ± 0.03
|
Suco
results <- summ_qf(dframe_qf_ufscar,variable_name = "Suco")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
48.76667
|
4.05010
|
48.77 ± 4.05
|
|
S
|
FD+H10
|
37.46667
|
3.35012
|
37.47 ± 3.35
|
|
S
|
H14+H10
|
45.63333
|
1.50444
|
45.63 ± 1.5
|
|
S
|
H152+H10
|
45.13333
|
2.58908
|
45.13 ± 2.59
|
|
S
|
H228+H10
|
39.86667
|
0.46188
|
39.87 ± 0.46
|
|
TRF
|
FD+BRS_PF
|
51.83333
|
2.10792
|
51.83 ± 2.11
|
|
TRF
|
FD+H10
|
51.20000
|
2.07846
|
51.2 ± 2.08
|
|
TRF
|
H14+H10
|
46.46667
|
1.76163
|
46.47 ± 1.76
|
|
TRF
|
H152+H10
|
45.00000
|
2.70000
|
45 ± 2.7
|
|
TRF
|
H228+H10
|
45.63333
|
3.13900
|
45.63 ± 3.14
|
|
TRV
|
FD+BRS_PF
|
45.13333
|
0.05774
|
45.13 ± 0.06
|
|
TRV
|
FD+H10
|
49.56667
|
0.75056
|
49.57 ± 0.75
|
|
TRV
|
H14+H10
|
38.63333
|
4.06735
|
38.63 ± 4.07
|
|
TRV
|
H152+H10
|
41.90000
|
1.60935
|
41.9 ± 1.61
|
|
TRV
|
H228+H10
|
43.43333
|
3.76475
|
43.43 ± 3.76
|
Acidez
results <- summ_qf(dframe_qf_ufscar,variable_name = "Acidez")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
5.29067
|
0.38493
|
5.29 ± 0.38
|
|
S
|
FD+H10
|
5.42190
|
0.28282
|
5.42 ± 0.28
|
|
S
|
H14+H10
|
5.01870
|
0.39804
|
5.02 ± 0.4
|
|
S
|
H152+H10
|
4.60677
|
0.51469
|
4.61 ± 0.51
|
|
S
|
H228+H10
|
5.06987
|
0.18971
|
5.07 ± 0.19
|
|
TRF
|
FD+BRS_PF
|
5.46163
|
0.26236
|
5.46 ± 0.26
|
|
TRF
|
FD+H10
|
5.14260
|
0.43562
|
5.14 ± 0.44
|
|
TRF
|
H14+H10
|
4.97363
|
0.08585
|
4.97 ± 0.09
|
|
TRF
|
H152+H10
|
5.16717
|
0.15245
|
5.17 ± 0.15
|
|
TRF
|
H228+H10
|
5.51613
|
0.34844
|
5.52 ± 0.35
|
|
TRV
|
FD+BRS_PF
|
5.15400
|
0.32922
|
5.15 ± 0.33
|
|
TRV
|
FD+H10
|
4.99920
|
0.20471
|
5 ± 0.2
|
|
TRV
|
H14+H10
|
5.00657
|
0.06020
|
5.01 ± 0.06
|
|
TRV
|
H152+H10
|
4.79327
|
0.15489
|
4.79 ± 0.15
|
|
TRV
|
H228+H10
|
5.50537
|
0.33448
|
5.51 ± 0.33
|
BrixCorr
results <- summ_qf(dframe_qf_ufscar,variable_name = "BrixCorr")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
8.90000
|
0.17321
|
8.9 ± 0.17
|
|
S
|
FD+H10
|
8.80000
|
0.10000
|
8.8 ± 0.1
|
|
S
|
H14+H10
|
8.83333
|
0.15275
|
8.83 ± 0.15
|
|
S
|
H152+H10
|
8.63333
|
0.05774
|
8.63 ± 0.06
|
|
S
|
H228+H10
|
9.56667
|
0.20817
|
9.57 ± 0.21
|
|
TRF
|
FD+BRS_PF
|
9.13333
|
0.15275
|
9.13 ± 0.15
|
|
TRF
|
FD+H10
|
9.16667
|
0.05774
|
9.17 ± 0.06
|
|
TRF
|
H14+H10
|
8.93333
|
0.11547
|
8.93 ± 0.12
|
|
TRF
|
H152+H10
|
8.90000
|
0.10000
|
8.9 ± 0.1
|
|
TRF
|
H228+H10
|
8.90000
|
0.00000
|
8.9 ± 0
|
|
TRV
|
FD+BRS_PF
|
8.83333
|
0.20817
|
8.83 ± 0.21
|
|
TRV
|
FD+H10
|
9.20000
|
0.10000
|
9.2 ± 0.1
|
|
TRV
|
H14+H10
|
8.86667
|
0.20817
|
8.87 ± 0.21
|
|
TRV
|
H152+H10
|
8.66667
|
0.05774
|
8.67 ± 0.06
|
|
TRV
|
H228+H10
|
8.96667
|
0.25166
|
8.97 ± 0.25
|
Ratio
results <- summ_qf(dframe_qf_ufscar,variable_name = "Ratio")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
1.70000
|
0.10000
|
1.7 ± 0.1
|
|
S
|
FD+H10
|
1.60000
|
0.10000
|
1.6 ± 0.1
|
|
S
|
H14+H10
|
1.73333
|
0.11547
|
1.73 ± 0.12
|
|
S
|
H152+H10
|
1.90000
|
0.20000
|
1.9 ± 0.2
|
|
S
|
H228+H10
|
1.90000
|
0.10000
|
1.9 ± 0.1
|
|
TRF
|
FD+BRS_PF
|
1.66667
|
0.05774
|
1.67 ± 0.06
|
|
TRF
|
FD+H10
|
1.76667
|
0.15275
|
1.77 ± 0.15
|
|
TRF
|
H14+H10
|
1.80000
|
0.00000
|
1.8 ± 0
|
|
TRF
|
H152+H10
|
1.75000
|
0.05000
|
1.75 ± 0.05
|
|
TRF
|
H228+H10
|
1.60000
|
0.10000
|
1.6 ± 0.1
|
|
TRV
|
FD+BRS_PF
|
1.73333
|
0.15275
|
1.73 ± 0.15
|
|
TRV
|
FD+H10
|
1.86667
|
0.05774
|
1.87 ± 0.06
|
|
TRV
|
H14+H10
|
1.76667
|
0.05774
|
1.77 ± 0.06
|
|
TRV
|
H152+H10
|
1.80000
|
0.10000
|
1.8 ± 0.1
|
|
TRV
|
H228+H10
|
1.63333
|
0.05774
|
1.63 ± 0.06
|
IT
results <- summ_qf(dframe_qf_ufscar,variable_name = "IT")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
1.77000
|
0.18358
|
1.77 ± 0.18
|
|
S
|
FD+H10
|
1.34333
|
0.11015
|
1.34 ± 0.11
|
|
S
|
H14+H10
|
1.64667
|
0.03512
|
1.65 ± 0.04
|
|
S
|
H152+H10
|
1.58667
|
0.09713
|
1.59 ± 0.1
|
|
S
|
H228+H10
|
1.55333
|
0.02309
|
1.55 ± 0.02
|
|
TRF
|
FD+BRS_PF
|
1.93000
|
0.05292
|
1.93 ± 0.05
|
|
TRF
|
FD+H10
|
1.91667
|
0.05508
|
1.92 ± 0.06
|
|
TRF
|
H14+H10
|
1.69000
|
0.08185
|
1.69 ± 0.08
|
|
TRF
|
H152+H10
|
1.63667
|
0.08505
|
1.64 ± 0.09
|
|
TRF
|
H228+H10
|
1.66000
|
0.11533
|
1.66 ± 0.12
|
|
TRV
|
FD+BRS_PF
|
1.62667
|
0.04509
|
1.63 ± 0.05
|
|
TRV
|
FD+H10
|
1.86000
|
0.02000
|
1.86 ± 0.02
|
|
TRV
|
H14+H10
|
1.39667
|
0.18148
|
1.4 ± 0.18
|
|
TRV
|
H152+H10
|
1.48000
|
0.06083
|
1.48 ± 0.06
|
|
TRV
|
H228+H10
|
1.60000
|
0.17349
|
1.6 ± 0.17
|
VitC
results <- summ_qf(dframe_qf_ufscar,variable_name = "VitC")
Gráfico de Barras
results$barp

results$barp_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
RegHidrico
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
S
|
FD+BRS_PF
|
28.54000
|
0.60622
|
28.54 ± 0.61
|
|
S
|
FD+H10
|
35.68000
|
1.14582
|
35.68 ± 1.15
|
|
S
|
H14+H10
|
32.45667
|
3.28649
|
32.46 ± 3.29
|
|
S
|
H152+H10
|
37.78333
|
1.63102
|
37.78 ± 1.63
|
|
S
|
H228+H10
|
33.36000
|
1.15763
|
33.36 ± 1.16
|
|
TRF
|
FD+BRS_PF
|
33.21000
|
1.82371
|
33.21 ± 1.82
|
|
TRF
|
FD+H10
|
32.76000
|
1.84984
|
32.76 ± 1.85
|
|
TRF
|
H14+H10
|
32.50667
|
1.82870
|
32.51 ± 1.83
|
|
TRF
|
H152+H10
|
34.24333
|
0.89131
|
34.24 ± 0.89
|
|
TRF
|
H228+H10
|
34.38333
|
1.21187
|
34.38 ± 1.21
|
|
TRV
|
FD+BRS_PF
|
33.45333
|
1.59080
|
33.45 ± 1.59
|
|
TRV
|
FD+H10
|
35.30333
|
0.89891
|
35.3 ± 0.9
|
|
TRV
|
H14+H10
|
38.17667
|
1.34500
|
38.18 ± 1.35
|
|
TRV
|
H152+H10
|
37.89000
|
1.84789
|
37.89 ± 1.85
|
|
TRV
|
H228+H10
|
35.40333
|
2.43097
|
35.4 ± 2.43
|
CCSM
Massa
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Massa")
models <- fit_qf_ccsm(dframe_qf_ccsm,'Massa')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
Manejo
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
Convencional
|
BRF_PF+FD
|
94.33333
|
11.50362
|
94.33 ± 11.5
|
|
Convencional
|
BRF_PF+H152
|
89.33333
|
10.01665
|
89.33 ± 10.02
|
|
Convencional
|
IAC10+FD
|
90.33333
|
4.72582
|
90.33 ± 4.73
|
|
Convencional
|
IAC10+H152
|
86.66667
|
5.13160
|
86.67 ± 5.13
|
|
Convencional
|
IAC5+FD
|
88.00000
|
12.12436
|
88 ± 12.12
|
|
Convencional
|
IAC5+H152
|
81.66667
|
2.30940
|
81.67 ± 2.31
|
|
Ecologico
|
BRF_PF+FD
|
77.33333
|
3.05505
|
77.33 ± 3.06
|
|
Ecologico
|
BRF_PF+H152
|
83.66667
|
0.57735
|
83.67 ± 0.58
|
|
Ecologico
|
IAC10+FD
|
95.00000
|
11.13553
|
95 ± 11.14
|
|
Ecologico
|
IAC10+H152
|
90.33333
|
2.51661
|
90.33 ± 2.52
|
|
Ecologico
|
IAC5+FD
|
80.33333
|
3.21455
|
80.33 ± 3.21
|
|
Ecologico
|
IAC5+H152
|
88.00000
|
3.46410
|
88 ± 3.46
|
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
|
Manejo
|
PortaEnxerto
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
FD
|
IAC10
|
95.00000
|
a
|
|
Convencional
|
FD
|
BRF_PF
|
94.33333
|
a
|
|
Convencional
|
FD
|
IAC10
|
90.33333
|
a
|
|
Ecologico
|
H152
|
IAC10
|
90.33333
|
a
|
|
Convencional
|
H152
|
BRF_PF
|
89.33333
|
a
|
|
Convencional
|
FD
|
IAC5
|
88.00000
|
a
|
|
Ecologico
|
H152
|
IAC5
|
88.00000
|
a
|
|
Convencional
|
H152
|
IAC10
|
86.66667
|
a
|
|
Ecologico
|
H152
|
BRF_PF
|
83.66667
|
a
|
|
Convencional
|
H152
|
IAC5
|
81.66667
|
a
|
|
Ecologico
|
FD
|
IAC5
|
80.33333
|
a
|
|
Ecologico
|
FD
|
BRF_PF
|
77.33333
|
a
|
summary(models$models$fit_full_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 61.4 61.36 1.235 0.2774
## PortaEnxerto 1 8.0 8.03 0.162 0.6912
## Copa 2 237.2 118.58 2.388 0.1133
## Manejo:PortaEnxerto 1 148.0 148.03 2.980 0.0971 .
## Manejo:Copa 2 377.4 188.69 3.799 0.0369 *
## PortaEnxerto:Copa 2 46.7 23.36 0.470 0.6304
## Manejo:PortaEnxerto:Copa 2 96.1 48.03 0.967 0.3945
## Residuals 24 1192.0 49.67
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
94.33333
|
4.06885
|
23.18426
|
0.00000
|
|
ManejoEcologico
|
-17.00000
|
5.75423
|
-2.95435
|
0.00691
|
|
PortaEnxertoH152
|
-5.00000
|
5.75423
|
-0.86893
|
0.39349
|
|
CopaIAC10
|
-4.00000
|
5.75423
|
-0.69514
|
0.49364
|
|
CopaIAC5
|
-6.33333
|
5.75423
|
-1.10064
|
0.28197
|
|
ManejoEcologico:PortaEnxertoH152
|
11.33333
|
8.13770
|
1.39269
|
0.17648
|
|
ManejoEcologico:CopaIAC10
|
21.66667
|
8.13770
|
2.66250
|
0.01363
|
|
ManejoEcologico:CopaIAC5
|
9.33333
|
8.13770
|
1.14692
|
0.26272
|
|
PortaEnxertoH152:CopaIAC10
|
1.33333
|
8.13770
|
0.16385
|
0.87122
|
|
PortaEnxertoH152:CopaIAC5
|
-1.33333
|
8.13770
|
-0.16385
|
0.87122
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC10
|
-12.33333
|
11.50845
|
-1.07168
|
0.29452
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC5
|
2.66667
|
11.50845
|
0.23171
|
0.81872
|
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
IAC10
|
92.66667
|
a
|
|
Convencional
|
BRF_PF
|
91.83333
|
a
|
|
Convencional
|
IAC10
|
88.50000
|
a
|
|
Convencional
|
IAC5
|
84.83333
|
a
|
|
Ecologico
|
IAC5
|
84.16667
|
a
|
|
Ecologico
|
BRF_PF
|
80.50000
|
a
|
[[2]]
|
Manejo
|
PortaEnxerto
|
Mean
|
Group
|
|
Convencional
|
FD
|
90.88889
|
a
|
|
Ecologico
|
H152
|
87.33333
|
a
|
|
Convencional
|
H152
|
85.88889
|
a
|
|
Ecologico
|
FD
|
84.22222
|
a
|
[[3]]
|
Copa
|
PortaEnxerto
|
Mean
|
Group
|
|
IAC10
|
FD
|
92.66667
|
a
|
|
IAC10
|
H152
|
88.50000
|
a
|
|
BRF_PF
|
H152
|
86.50000
|
a
|
|
BRF_PF
|
FD
|
85.83333
|
a
|
|
IAC5
|
H152
|
84.83333
|
a
|
|
IAC5
|
FD
|
84.16667
|
a
|
summary(models$models$fit_noInt3_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 61.4 61.36 1.239 0.2759
## Copa 2 237.2 118.58 2.394 0.1111
## PortaEnxerto 1 8.0 8.03 0.162 0.6906
## Manejo:Copa 2 377.4 188.69 3.809 0.0354 *
## Manejo:PortaEnxerto 1 148.0 148.03 2.988 0.0957 .
## Copa:PortaEnxerto 2 46.7 23.36 0.472 0.6293
## Residuals 26 1288.1 49.54
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
93.52778
|
3.70962
|
25.21223
|
0.00000
|
|
ManejoEcologico
|
-15.38889
|
4.69234
|
-3.27958
|
0.00296
|
|
CopaIAC10
|
-0.91667
|
4.97698
|
-0.18418
|
0.85530
|
|
CopaIAC5
|
-7.00000
|
4.97698
|
-1.40648
|
0.17142
|
|
PortaEnxertoH152
|
-3.38889
|
4.69234
|
-0.72222
|
0.47661
|
|
ManejoEcologico:CopaIAC10
|
15.50000
|
5.74692
|
2.69710
|
0.01211
|
|
ManejoEcologico:CopaIAC5
|
10.66667
|
5.74692
|
1.85607
|
0.07481
|
|
ManejoEcologico:PortaEnxertoH152
|
8.11111
|
4.69234
|
1.72859
|
0.09574
|
|
CopaIAC10:PortaEnxertoH152
|
-4.83333
|
5.74692
|
-0.84103
|
0.40800
|
|
CopaIAC5:PortaEnxertoH152
|
0.00000
|
5.74692
|
0.00000
|
1.00000
|
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
IAC10
|
92.66667
|
a
|
|
Convencional
|
BRF_PF
|
91.83333
|
a
|
|
Convencional
|
IAC10
|
88.50000
|
a
|
|
Convencional
|
IAC5
|
84.83333
|
a
|
|
Ecologico
|
IAC5
|
84.16667
|
a
|
|
Ecologico
|
BRF_PF
|
80.50000
|
a
|
[[2]]
|
Manejo
|
Mean
|
Group
|
|
Convencional
|
88.38889
|
a
|
|
Ecologico
|
85.77778
|
a
|
[[3]]
|
Copa
|
Mean
|
Group
|
|
IAC10
|
90.58333
|
a
|
|
BRF_PF
|
86.16667
|
a
|
|
IAC5
|
84.50000
|
a
|
[[4]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
87.55556
|
a
|
|
H152
|
86.61111
|
a
|
summary(models$models$fit_CopPort_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 61.4 61.36 0.981 0.330
## PortaEnxerto 1 8.0 8.03 0.128 0.723
## Copa 2 237.2 118.58 1.896 0.168
## PortaEnxerto:Copa 2 46.7 23.36 0.374 0.692
## Residuals 29 1813.5 62.53
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
87.13889
|
3.48702
|
24.98951
|
0.00000
|
|
ManejoEcologico
|
-2.61111
|
2.63594
|
-0.99058
|
0.33008
|
|
PortaEnxertoH152
|
0.66667
|
4.56558
|
0.14602
|
0.88492
|
|
CopaIAC10
|
6.83333
|
4.56558
|
1.49671
|
0.14527
|
|
CopaIAC5
|
-1.66667
|
4.56558
|
-0.36505
|
0.71772
|
|
PortaEnxertoH152:CopaIAC10
|
-4.83333
|
6.45670
|
-0.74858
|
0.46014
|
|
PortaEnxertoH152:CopaIAC5
|
0.00000
|
6.45670
|
0.00000
|
1.00000
|
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
|
Manejo
|
Mean
|
Group
|
|
Convencional
|
88.38889
|
a
|
|
Ecologico
|
85.77778
|
a
|
[[2]]
|
Copa
|
Mean
|
Group
|
|
IAC10
|
90.58333
|
a
|
|
BRF_PF
|
86.16667
|
a
|
|
IAC5
|
84.50000
|
a
|
[[3]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
87.55556
|
a
|
|
H152
|
86.61111
|
a
|
summary(models$models$fit_mainEffect_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 61.4 61.36 1.023 0.320
## PortaEnxerto 1 8.0 8.03 0.134 0.717
## Copa 2 237.2 118.58 1.976 0.156
## Residuals 31 1860.2 60.01
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
87.94444
|
2.88690
|
30.46326
|
0.00000
|
|
ManejoEcologico
|
-2.61111
|
2.58212
|
-1.01123
|
0.31974
|
|
PortaEnxertoH152
|
-0.94444
|
2.58212
|
-0.36576
|
0.71703
|
|
CopaIAC10
|
4.41667
|
3.16244
|
1.39660
|
0.17245
|
|
CopaIAC5
|
-1.66667
|
3.16244
|
-0.52702
|
0.60193
|
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))
Altura
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Altura")
models <- fit_qf_ccsm(dframe_qf_ccsm,'Altura')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
Manejo
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
Convencional
|
BRF_PF+FD
|
6.10000
|
0.40000
|
6.1 ± 0.4
|
|
Convencional
|
BRF_PF+H152
|
5.96667
|
0.25166
|
5.97 ± 0.25
|
|
Convencional
|
IAC10+FD
|
6.06667
|
0.15275
|
6.07 ± 0.15
|
|
Convencional
|
IAC10+H152
|
5.56667
|
0.11547
|
5.57 ± 0.12
|
|
Convencional
|
IAC5+FD
|
5.90000
|
0.36056
|
5.9 ± 0.36
|
|
Convencional
|
IAC5+H152
|
5.80000
|
0.10000
|
5.8 ± 0.1
|
|
Ecologico
|
BRF_PF+FD
|
5.70000
|
0.10000
|
5.7 ± 0.1
|
|
Ecologico
|
BRF_PF+H152
|
5.83333
|
0.05774
|
5.83 ± 0.06
|
|
Ecologico
|
IAC10+FD
|
6.13333
|
0.15275
|
6.13 ± 0.15
|
|
Ecologico
|
IAC10+H152
|
6.03333
|
0.11547
|
6.03 ± 0.12
|
|
Ecologico
|
IAC5+FD
|
5.90000
|
0.10000
|
5.9 ± 0.1
|
|
Ecologico
|
IAC5+H152
|
6.03333
|
0.05774
|
6.03 ± 0.06
|
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
|
Manejo
|
PortaEnxerto
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
FD
|
IAC10
|
6.13333
|
a
|
|
Convencional
|
FD
|
BRF_PF
|
6.10000
|
a
|
|
Convencional
|
FD
|
IAC10
|
6.06667
|
a
|
|
Ecologico
|
H152
|
IAC10
|
6.03333
|
a
|
|
Ecologico
|
H152
|
IAC5
|
6.03333
|
a
|
|
Convencional
|
H152
|
BRF_PF
|
5.96667
|
a
|
|
Ecologico
|
FD
|
IAC5
|
5.90000
|
a
|
|
Convencional
|
FD
|
IAC5
|
5.90000
|
a
|
|
Ecologico
|
H152
|
BRF_PF
|
5.83333
|
a
|
|
Convencional
|
H152
|
IAC5
|
5.80000
|
a
|
|
Ecologico
|
FD
|
BRF_PF
|
5.70000
|
a
|
|
Convencional
|
H152
|
IAC10
|
5.56667
|
a
|
summary(models$models$fit_full_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0136 0.01361 0.353 0.55825
## PortaEnxerto 1 0.0803 0.08028 2.079 0.16225
## Copa 2 0.0172 0.00861 0.223 0.80174
## Manejo:PortaEnxerto 1 0.2025 0.20250 5.245 0.03110 *
## Manejo:Copa 2 0.4539 0.22694 5.878 0.00836 **
## PortaEnxerto:Copa 2 0.1906 0.09528 2.468 0.10602
## Manejo:PortaEnxerto:Copa 2 0.0117 0.00583 0.151 0.86059
## Residuals 24 0.9267 0.03861
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
6.10000
|
0.11345
|
53.76929
|
0.00000
|
|
ManejoEcologico
|
-0.40000
|
0.16044
|
-2.49316
|
0.01996
|
|
PortaEnxertoH152
|
-0.13333
|
0.16044
|
-0.83105
|
0.41414
|
|
CopaIAC10
|
-0.03333
|
0.16044
|
-0.20776
|
0.83717
|
|
CopaIAC5
|
-0.20000
|
0.16044
|
-1.24658
|
0.22458
|
|
ManejoEcologico:PortaEnxertoH152
|
0.26667
|
0.22690
|
1.17529
|
0.25141
|
|
ManejoEcologico:CopaIAC10
|
0.46667
|
0.22690
|
2.05675
|
0.05074
|
|
ManejoEcologico:CopaIAC5
|
0.40000
|
0.22690
|
1.76293
|
0.09064
|
|
PortaEnxertoH152:CopaIAC10
|
-0.36667
|
0.22690
|
-1.61602
|
0.11916
|
|
PortaEnxertoH152:CopaIAC5
|
0.03333
|
0.22690
|
0.14691
|
0.88443
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC10
|
0.13333
|
0.32088
|
0.41553
|
0.68145
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC5
|
-0.03333
|
0.32088
|
-0.10388
|
0.91813
|
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
IAC10
|
6.08333
|
a
|
|
Convencional
|
BRF_PF
|
6.03333
|
a
|
|
Ecologico
|
IAC5
|
5.96667
|
a
|
|
Convencional
|
IAC5
|
5.85000
|
a
|
|
Convencional
|
IAC10
|
5.81667
|
a
|
|
Ecologico
|
BRF_PF
|
5.76667
|
a
|
[[2]]
|
Manejo
|
PortaEnxerto
|
Mean
|
Group
|
|
Convencional
|
FD
|
6.02222
|
a
|
|
Ecologico
|
H152
|
5.96667
|
a
|
|
Ecologico
|
FD
|
5.91111
|
a
|
|
Convencional
|
H152
|
5.77778
|
a
|
[[3]]
|
Copa
|
PortaEnxerto
|
Mean
|
Group
|
|
IAC10
|
FD
|
6.10000
|
a
|
|
IAC5
|
H152
|
5.91667
|
a
|
|
BRF_PF
|
FD
|
5.90000
|
a
|
|
BRF_PF
|
H152
|
5.90000
|
a
|
|
IAC5
|
FD
|
5.90000
|
a
|
|
IAC10
|
H152
|
5.80000
|
a
|
summary(models$models$fit_noInt3_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0136 0.01361 0.377 0.54447
## Copa 2 0.0172 0.00861 0.239 0.78943
## PortaEnxerto 1 0.0803 0.08028 2.224 0.14788
## Manejo:Copa 2 0.4539 0.22694 6.288 0.00592 **
## Manejo:PortaEnxerto 1 0.2025 0.20250 5.611 0.02556 *
## Copa:PortaEnxerto 2 0.1906 0.09528 2.640 0.09040 .
## Residuals 26 0.9383 0.03609
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
6.10833
|
0.10012
|
61.00734
|
0.00000
|
|
ManejoEcologico
|
-0.41667
|
0.12665
|
-3.28994
|
0.00288
|
|
CopaIAC10
|
-0.06667
|
0.13433
|
-0.49629
|
0.62386
|
|
CopaIAC5
|
-0.19167
|
0.13433
|
-1.42682
|
0.16553
|
|
PortaEnxertoH152
|
-0.15000
|
0.12665
|
-1.18438
|
0.24698
|
|
ManejoEcologico:CopaIAC10
|
0.53333
|
0.15511
|
3.43837
|
0.00198
|
|
ManejoEcologico:CopaIAC5
|
0.38333
|
0.15511
|
2.47133
|
0.02033
|
|
ManejoEcologico:PortaEnxertoH152
|
0.30000
|
0.12665
|
2.36876
|
0.02556
|
|
CopaIAC10:PortaEnxertoH152
|
-0.30000
|
0.15511
|
-1.93408
|
0.06406
|
|
CopaIAC5:PortaEnxertoH152
|
0.01667
|
0.15511
|
0.10745
|
0.91526
|
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
IAC10
|
6.08333
|
a
|
|
Convencional
|
BRF_PF
|
6.03333
|
a
|
|
Ecologico
|
IAC5
|
5.96667
|
a
|
|
Convencional
|
IAC5
|
5.85000
|
a
|
|
Convencional
|
IAC10
|
5.81667
|
a
|
|
Ecologico
|
BRF_PF
|
5.76667
|
a
|
[[2]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
5.93889
|
a
|
|
Convencional
|
5.90000
|
a
|
[[3]]
|
Copa
|
Mean
|
Group
|
|
IAC10
|
5.95000
|
a
|
|
IAC5
|
5.90833
|
a
|
|
BRF_PF
|
5.90000
|
a
|
[[4]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
5.96667
|
a
|
|
H152
|
5.87222
|
a
|
summary(models$models$fit_CopPort_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0136 0.01361 0.248 0.623
## PortaEnxerto 1 0.0803 0.08028 1.460 0.237
## Copa 2 0.0172 0.00861 0.157 0.856
## PortaEnxerto:Copa 2 0.1906 0.09528 1.733 0.195
## Residuals 29 1.5947 0.05499
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
5.88056
|
0.10340
|
56.86919
|
0.00000
|
|
ManejoEcologico
|
0.03889
|
0.07817
|
0.49751
|
0.62258
|
|
PortaEnxertoH152
|
0.00000
|
0.13539
|
0.00000
|
1.00000
|
|
CopaIAC10
|
0.20000
|
0.13539
|
1.47723
|
0.15039
|
|
CopaIAC5
|
0.00000
|
0.13539
|
0.00000
|
1.00000
|
|
PortaEnxertoH152:CopaIAC10
|
-0.30000
|
0.19147
|
-1.56684
|
0.12800
|
|
PortaEnxertoH152:CopaIAC5
|
0.01667
|
0.19147
|
0.08705
|
0.93123
|
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
5.93889
|
a
|
|
Convencional
|
5.90000
|
a
|
[[2]]
|
Copa
|
Mean
|
Group
|
|
IAC10
|
5.95000
|
a
|
|
IAC5
|
5.90833
|
a
|
|
BRF_PF
|
5.90000
|
a
|
[[3]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
5.96667
|
a
|
|
H152
|
5.87222
|
a
|
summary(models$models$fit_mainEffect_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0136 0.01361 0.236 0.630
## PortaEnxerto 1 0.0803 0.08028 1.394 0.247
## Copa 2 0.0172 0.00861 0.150 0.862
## Residuals 31 1.7853 0.05759
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
5.92778
|
0.08943
|
66.28055
|
0.00000
|
|
ManejoEcologico
|
0.03889
|
0.07999
|
0.48615
|
0.63028
|
|
PortaEnxertoH152
|
-0.09444
|
0.07999
|
-1.18066
|
0.24672
|
|
CopaIAC10
|
0.05000
|
0.09797
|
0.51036
|
0.61342
|
|
CopaIAC5
|
0.00833
|
0.09797
|
0.08506
|
0.93276
|
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))
AL
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "AL")
models <- fit_qf_ccsm(dframe_qf_ccsm,'AL')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
Manejo
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
Convencional
|
BRF_PF+FD
|
1.09000
|
0.01732
|
1.09 ± 0.02
|
|
Convencional
|
BRF_PF+H152
|
1.10667
|
0.02082
|
1.11 ± 0.02
|
|
Convencional
|
IAC10+FD
|
1.12333
|
0.01155
|
1.12 ± 0.01
|
|
Convencional
|
IAC10+H152
|
1.09333
|
0.01155
|
1.09 ± 0.01
|
|
Convencional
|
IAC5+FD
|
1.08667
|
0.02082
|
1.09 ± 0.02
|
|
Convencional
|
IAC5+H152
|
1.10000
|
0.01000
|
1.1 ± 0.01
|
|
Ecologico
|
BRF_PF+FD
|
1.10667
|
0.03055
|
1.11 ± 0.03
|
|
Ecologico
|
BRF_PF+H152
|
1.10667
|
0.01528
|
1.11 ± 0.02
|
|
Ecologico
|
IAC10+FD
|
1.10333
|
0.01155
|
1.1 ± 0.01
|
|
Ecologico
|
IAC10+H152
|
1.15333
|
0.03786
|
1.15 ± 0.04
|
|
Ecologico
|
IAC5+FD
|
1.11333
|
0.02082
|
1.11 ± 0.02
|
|
Ecologico
|
IAC5+H152
|
1.11000
|
0.02000
|
1.11 ± 0.02
|
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
|
Manejo
|
PortaEnxerto
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
H152
|
IAC10
|
1.15333
|
a
|
|
Convencional
|
FD
|
IAC10
|
1.12333
|
ab
|
|
Ecologico
|
FD
|
IAC5
|
1.11333
|
ab
|
|
Ecologico
|
H152
|
IAC5
|
1.11000
|
ab
|
|
Convencional
|
H152
|
BRF_PF
|
1.10667
|
ab
|
|
Ecologico
|
FD
|
BRF_PF
|
1.10667
|
ab
|
|
Ecologico
|
H152
|
BRF_PF
|
1.10667
|
ab
|
|
Ecologico
|
FD
|
IAC10
|
1.10333
|
ab
|
|
Convencional
|
H152
|
IAC5
|
1.10000
|
ab
|
|
Convencional
|
H152
|
IAC10
|
1.09333
|
ab
|
|
Convencional
|
FD
|
BRF_PF
|
1.09000
|
b
|
|
Convencional
|
FD
|
IAC5
|
1.08667
|
b
|
summary(models$models$fit_full_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.002178 0.0021778 5.124 0.0329 *
## PortaEnxerto 1 0.000544 0.0005444 1.281 0.2689
## Copa 2 0.002006 0.0010028 2.359 0.1160
## Manejo:PortaEnxerto 1 0.000544 0.0005444 1.281 0.2689
## Manejo:Copa 2 0.000239 0.0001194 0.281 0.7574
## PortaEnxerto:Copa 2 0.000039 0.0000194 0.046 0.9554
## Manejo:PortaEnxerto:Copa 2 0.004672 0.0023361 5.497 0.0108 *
## Residuals 24 0.010200 0.0004250
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.09000
|
0.01190
|
91.57832
|
0.00000
|
|
ManejoEcologico
|
0.01667
|
0.01683
|
0.99015
|
0.33198
|
|
PortaEnxertoH152
|
0.01667
|
0.01683
|
0.99015
|
0.33198
|
|
CopaIAC10
|
0.03333
|
0.01683
|
1.98030
|
0.05925
|
|
CopaIAC5
|
-0.00333
|
0.01683
|
-0.19803
|
0.84469
|
|
ManejoEcologico:PortaEnxertoH152
|
-0.01667
|
0.02380
|
-0.70014
|
0.49057
|
|
ManejoEcologico:CopaIAC10
|
-0.03667
|
0.02380
|
-1.54031
|
0.13657
|
|
ManejoEcologico:CopaIAC5
|
0.01000
|
0.02380
|
0.42008
|
0.67816
|
|
PortaEnxertoH152:CopaIAC10
|
-0.04667
|
0.02380
|
-1.96039
|
0.06166
|
|
PortaEnxertoH152:CopaIAC5
|
-0.00333
|
0.02380
|
-0.14003
|
0.88981
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC10
|
0.09667
|
0.03367
|
2.87143
|
0.00841
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC5
|
0.00000
|
0.03367
|
0.00000
|
1.00000
|
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
IAC10
|
1.12833
|
a
|
|
Ecologico
|
IAC5
|
1.11167
|
a
|
|
Convencional
|
IAC10
|
1.10833
|
a
|
|
Ecologico
|
BRF_PF
|
1.10667
|
a
|
|
Convencional
|
BRF_PF
|
1.09833
|
a
|
|
Convencional
|
IAC5
|
1.09333
|
a
|
[[2]]
|
Manejo
|
PortaEnxerto
|
Mean
|
Group
|
|
Ecologico
|
H152
|
1.12333
|
a
|
|
Ecologico
|
FD
|
1.10778
|
a
|
|
Convencional
|
FD
|
1.10000
|
a
|
|
Convencional
|
H152
|
1.10000
|
a
|
[[3]]
|
Copa
|
PortaEnxerto
|
Mean
|
Group
|
|
IAC10
|
H152
|
1.12333
|
a
|
|
IAC10
|
FD
|
1.11333
|
a
|
|
BRF_PF
|
H152
|
1.10667
|
a
|
|
IAC5
|
H152
|
1.10500
|
a
|
|
IAC5
|
FD
|
1.10000
|
a
|
|
BRF_PF
|
FD
|
1.09833
|
a
|
summary(models$models$fit_noInt3_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.002178 0.0021778 3.807 0.0619 .
## Copa 2 0.002006 0.0010028 1.753 0.1931
## PortaEnxerto 1 0.000544 0.0005444 0.952 0.3382
## Manejo:Copa 2 0.000239 0.0001194 0.209 0.8129
## Manejo:PortaEnxerto 1 0.000544 0.0005444 0.952 0.3382
## Copa:PortaEnxerto 2 0.000039 0.0000194 0.034 0.9666
## Residuals 26 0.014872 0.0005720
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.09806
|
0.01261
|
87.11126
|
0.00000
|
|
ManejoEcologico
|
0.00056
|
0.01594
|
0.03484
|
0.97247
|
|
CopaIAC10
|
0.00917
|
0.01691
|
0.54203
|
0.59241
|
|
CopaIAC5
|
-0.00333
|
0.01691
|
-0.19710
|
0.84528
|
|
PortaEnxertoH152
|
0.00056
|
0.01594
|
0.03484
|
0.97247
|
|
ManejoEcologico:CopaIAC10
|
0.01167
|
0.01953
|
0.59744
|
0.55538
|
|
ManejoEcologico:CopaIAC5
|
0.01000
|
0.01953
|
0.51209
|
0.61291
|
|
ManejoEcologico:PortaEnxertoH152
|
0.01556
|
0.01594
|
0.97561
|
0.33825
|
|
CopaIAC10:PortaEnxertoH152
|
0.00167
|
0.01953
|
0.08535
|
0.93264
|
|
CopaIAC5:PortaEnxertoH152
|
-0.00333
|
0.01953
|
-0.17070
|
0.86578
|
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
IAC10
|
1.12833
|
a
|
|
Ecologico
|
IAC5
|
1.11167
|
a
|
|
Convencional
|
IAC10
|
1.10833
|
a
|
|
Ecologico
|
BRF_PF
|
1.10667
|
a
|
|
Convencional
|
BRF_PF
|
1.09833
|
a
|
|
Convencional
|
IAC5
|
1.09333
|
a
|
[[2]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
1.11556
|
a
|
|
Convencional
|
1.10000
|
a
|
[[3]]
|
Copa
|
Mean
|
Group
|
|
IAC10
|
1.11833
|
a
|
|
BRF_PF
|
1.10250
|
a
|
|
IAC5
|
1.10250
|
a
|
[[4]]
|
PortaEnxerto
|
Mean
|
Group
|
|
H152
|
1.11167
|
a
|
|
FD
|
1.10389
|
a
|
summary(models$models$fit_CopPort_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.002178 0.0021778 4.034 0.054 .
## PortaEnxerto 1 0.000544 0.0005444 1.009 0.324
## Copa 2 0.002006 0.0010028 1.858 0.174
## PortaEnxerto:Copa 2 0.000039 0.0000194 0.036 0.965
## Residuals 29 0.015656 0.0005398
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.09056
|
0.01025
|
106.44243
|
0.00000
|
|
ManejoEcologico
|
0.01556
|
0.00774
|
2.00850
|
0.05398
|
|
PortaEnxertoH152
|
0.00833
|
0.01341
|
0.62122
|
0.53931
|
|
CopaIAC10
|
0.01500
|
0.01341
|
1.11819
|
0.27266
|
|
CopaIAC5
|
0.00167
|
0.01341
|
0.12424
|
0.90198
|
|
PortaEnxertoH152:CopaIAC10
|
0.00167
|
0.01897
|
0.08785
|
0.93060
|
|
PortaEnxertoH152:CopaIAC5
|
-0.00333
|
0.01897
|
-0.17571
|
0.86175
|
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
1.11556
|
a
|
|
Convencional
|
1.10000
|
b
|
[[2]]
|
Copa
|
Mean
|
Group
|
|
IAC10
|
1.11833
|
a
|
|
BRF_PF
|
1.10250
|
a
|
|
IAC5
|
1.10250
|
a
|
[[3]]
|
PortaEnxerto
|
Mean
|
Group
|
|
H152
|
1.11167
|
a
|
|
FD
|
1.10389
|
a
|
summary(models$models$fit_mainEffect_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.002178 0.0021778 4.302 0.0465 *
## PortaEnxerto 1 0.000544 0.0005444 1.075 0.3078
## Copa 2 0.002006 0.0010028 1.981 0.1551
## Residuals 31 0.015694 0.0005063
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.09083
|
0.00839
|
130.08659
|
0.00000
|
|
ManejoEcologico
|
0.01556
|
0.00750
|
2.07403
|
0.04647
|
|
PortaEnxertoH152
|
0.00778
|
0.00750
|
1.03701
|
0.30775
|
|
CopaIAC10
|
0.01583
|
0.00919
|
1.72368
|
0.09473
|
|
CopaIAC5
|
0.00000
|
0.00919
|
0.00000
|
1.00000
|
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))
Acidez
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Acidez")
models <- fit_qf_ccsm(dframe_qf_ccsm,'Acidez')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
Manejo
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
Convencional
|
BRF_PF+FD
|
5.79277
|
0.21457
|
5.79 ± 0.21
|
|
Convencional
|
BRF_PF+H152
|
4.73127
|
1.12089
|
4.73 ± 1.12
|
|
Convencional
|
IAC10+FD
|
5.78470
|
0.14056
|
5.78 ± 0.14
|
|
Convencional
|
IAC10+H152
|
4.64587
|
0.14255
|
4.65 ± 0.14
|
|
Convencional
|
IAC5+FD
|
5.50670
|
0.55673
|
5.51 ± 0.56
|
|
Convencional
|
IAC5+H152
|
5.28863
|
0.23447
|
5.29 ± 0.23
|
|
Ecologico
|
BRF_PF+FD
|
5.43200
|
0.25991
|
5.43 ± 0.26
|
|
Ecologico
|
BRF_PF+H152
|
5.71203
|
0.41669
|
5.71 ± 0.42
|
|
Ecologico
|
IAC10+FD
|
5.43063
|
0.73447
|
5.43 ± 0.73
|
|
Ecologico
|
IAC10+H152
|
5.57807
|
0.10602
|
5.58 ± 0.11
|
|
Ecologico
|
IAC5+FD
|
5.59083
|
0.15735
|
5.59 ± 0.16
|
|
Ecologico
|
IAC5+H152
|
5.64403
|
0.08699
|
5.64 ± 0.09
|
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
|
Manejo
|
PortaEnxerto
|
Copa
|
Mean
|
Group
|
|
Convencional
|
FD
|
BRF_PF
|
5.79277
|
a
|
|
Convencional
|
FD
|
IAC10
|
5.78470
|
a
|
|
Ecologico
|
H152
|
BRF_PF
|
5.71203
|
a
|
|
Ecologico
|
H152
|
IAC5
|
5.64403
|
a
|
|
Ecologico
|
FD
|
IAC5
|
5.59083
|
a
|
|
Ecologico
|
H152
|
IAC10
|
5.57807
|
a
|
|
Convencional
|
FD
|
IAC5
|
5.50670
|
a
|
|
Ecologico
|
FD
|
BRF_PF
|
5.43200
|
a
|
|
Ecologico
|
FD
|
IAC10
|
5.43063
|
a
|
|
Convencional
|
H152
|
IAC5
|
5.28863
|
a
|
|
Convencional
|
H152
|
BRF_PF
|
4.73127
|
a
|
|
Convencional
|
H152
|
IAC10
|
4.64587
|
a
|
summary(models$models$fit_full_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.670 0.6705 3.178 0.08728 .
## PortaEnxerto 1 0.939 0.9387 4.449 0.04553 *
## Copa 2 0.133 0.0666 0.316 0.73232
## Manejo:PortaEnxerto 1 2.101 2.1011 9.960 0.00427 **
## Manejo:Copa 2 0.013 0.0067 0.032 0.96882
## PortaEnxerto:Copa 2 0.277 0.1384 0.656 0.52790
## Manejo:PortaEnxerto:Copa 2 0.545 0.2723 1.291 0.29344
## Residuals 24 5.063 0.2110
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
5.79277
|
0.26518
|
21.84430
|
0.00000
|
|
ManejoEcologico
|
-0.36077
|
0.37503
|
-0.96197
|
0.34566
|
|
PortaEnxertoH152
|
-1.06150
|
0.37503
|
-2.83046
|
0.00925
|
|
CopaIAC10
|
-0.00807
|
0.37503
|
-0.02151
|
0.98302
|
|
CopaIAC5
|
-0.28607
|
0.37503
|
-0.76279
|
0.45302
|
|
ManejoEcologico:PortaEnxertoH152
|
1.34153
|
0.53037
|
2.52944
|
0.01840
|
|
ManejoEcologico:CopaIAC10
|
0.00670
|
0.53037
|
0.01263
|
0.99003
|
|
ManejoEcologico:CopaIAC5
|
0.44490
|
0.53037
|
0.83885
|
0.40983
|
|
PortaEnxertoH152:CopaIAC10
|
-0.07733
|
0.53037
|
-0.14581
|
0.88529
|
|
PortaEnxertoH152:CopaIAC5
|
0.84343
|
0.53037
|
1.59028
|
0.12486
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC10
|
-0.05527
|
0.75005
|
-0.07368
|
0.94187
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC5
|
-1.07027
|
0.75005
|
-1.42692
|
0.16648
|
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
IAC5
|
5.61743
|
a
|
|
Ecologico
|
BRF_PF
|
5.57202
|
a
|
|
Ecologico
|
IAC10
|
5.50435
|
a
|
|
Convencional
|
IAC5
|
5.39767
|
a
|
|
Convencional
|
BRF_PF
|
5.26202
|
a
|
|
Convencional
|
IAC10
|
5.21528
|
a
|
[[2]]
|
Manejo
|
PortaEnxerto
|
Mean
|
Group
|
|
Convencional
|
FD
|
5.69472
|
a
|
|
Ecologico
|
H152
|
5.64471
|
a
|
|
Ecologico
|
FD
|
5.48449
|
ab
|
|
Convencional
|
H152
|
4.88859
|
b
|
[[3]]
|
Copa
|
PortaEnxerto
|
Mean
|
Group
|
|
BRF_PF
|
FD
|
5.61238
|
a
|
|
IAC10
|
FD
|
5.60767
|
a
|
|
IAC5
|
FD
|
5.54877
|
a
|
|
IAC5
|
H152
|
5.46633
|
a
|
|
BRF_PF
|
H152
|
5.22165
|
a
|
|
IAC10
|
H152
|
5.11197
|
a
|
summary(models$models$fit_noInt3_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.670 0.6705 3.109 0.08963 .
## Copa 2 0.133 0.0666 0.309 0.73705
## PortaEnxerto 1 0.939 0.9387 4.352 0.04692 *
## Manejo:Copa 2 0.013 0.0067 0.031 0.96949
## Manejo:PortaEnxerto 1 2.101 2.1011 9.742 0.00438 **
## Copa:PortaEnxerto 2 0.277 0.1384 0.642 0.53448
## Residuals 26 5.608 0.2157
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
5.69897
|
0.24477
|
23.28270
|
0.00000
|
|
ManejoEcologico
|
-0.17318
|
0.30962
|
-0.55933
|
0.58072
|
|
CopaIAC10
|
0.00575
|
0.32840
|
0.01751
|
0.98616
|
|
CopaIAC5
|
-0.01850
|
0.32840
|
-0.05633
|
0.95551
|
|
PortaEnxertoH152
|
-0.87391
|
0.30962
|
-2.82257
|
0.00901
|
|
ManejoEcologico:CopaIAC10
|
-0.02093
|
0.37920
|
-0.05520
|
0.95640
|
|
ManejoEcologico:CopaIAC5
|
-0.09023
|
0.37920
|
-0.23796
|
0.81378
|
|
ManejoEcologico:PortaEnxertoH152
|
0.96636
|
0.30962
|
3.12114
|
0.00438
|
|
CopaIAC10:PortaEnxertoH152
|
-0.10497
|
0.37920
|
-0.27681
|
0.78411
|
|
CopaIAC5:PortaEnxertoH152
|
0.30830
|
0.37920
|
0.81303
|
0.42359
|
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
IAC5
|
5.61743
|
a
|
|
Ecologico
|
BRF_PF
|
5.57202
|
a
|
|
Ecologico
|
IAC10
|
5.50435
|
a
|
|
Convencional
|
IAC5
|
5.39767
|
a
|
|
Convencional
|
BRF_PF
|
5.26202
|
a
|
|
Convencional
|
IAC10
|
5.21528
|
a
|
[[2]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
5.56460
|
a
|
|
Convencional
|
5.29166
|
a
|
[[3]]
|
Copa
|
Mean
|
Group
|
|
IAC5
|
5.50755
|
a
|
|
BRF_PF
|
5.41702
|
a
|
|
IAC10
|
5.35982
|
a
|
[[4]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
5.58961
|
a
|
|
H152
|
5.26665
|
a
|
summary(models$models$fit_CopPort_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.670 0.6705 2.518 0.1234
## PortaEnxerto 1 0.939 0.9387 3.525 0.0705 .
## Copa 2 0.133 0.0666 0.250 0.7804
## PortaEnxerto:Copa 2 0.277 0.1384 0.520 0.6001
## Residuals 29 7.722 0.2663
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
5.47591
|
0.22755
|
24.06467
|
0.00000
|
|
ManejoEcologico
|
0.27294
|
0.17201
|
1.58678
|
0.12341
|
|
PortaEnxertoH152
|
-0.39073
|
0.29793
|
-1.31148
|
0.19998
|
|
CopaIAC10
|
-0.00472
|
0.29793
|
-0.01583
|
0.98748
|
|
CopaIAC5
|
-0.06362
|
0.29793
|
-0.21353
|
0.83241
|
|
PortaEnxertoH152:CopaIAC10
|
-0.10497
|
0.42134
|
-0.24913
|
0.80502
|
|
PortaEnxertoH152:CopaIAC5
|
0.30830
|
0.42134
|
0.73171
|
0.47022
|
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
5.56460
|
a
|
|
Convencional
|
5.29166
|
a
|
[[2]]
|
Copa
|
Mean
|
Group
|
|
IAC5
|
5.50755
|
a
|
|
BRF_PF
|
5.41702
|
a
|
|
IAC10
|
5.35982
|
a
|
[[3]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
5.58961
|
a
|
|
H152
|
5.26665
|
a
|
summary(models$models$fit_mainEffect_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.670 0.6705 2.598 0.1171
## PortaEnxerto 1 0.939 0.9387 3.638 0.0658 .
## Copa 2 0.133 0.0666 0.258 0.7742
## Residuals 31 7.999 0.2580
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
5.44202
|
0.18931
|
28.74624
|
0.00000
|
|
ManejoEcologico
|
0.27294
|
0.16933
|
1.61194
|
0.11711
|
|
PortaEnxertoH152
|
-0.32296
|
0.16933
|
-1.90730
|
0.06578
|
|
CopaIAC10
|
-0.05720
|
0.20738
|
-0.27582
|
0.78452
|
|
CopaIAC5
|
0.09053
|
0.20738
|
0.43655
|
0.66546
|
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))
BrixCorr
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "BrixCorr")
models <- fit_qf_ccsm(dframe_qf_ccsm,'BrixCorr')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
Manejo
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
Convencional
|
BRF_PF+FD
|
9.03333
|
0.15275
|
9.03 ± 0.15
|
|
Convencional
|
BRF_PF+H152
|
9.03333
|
0.15275
|
9.03 ± 0.15
|
|
Convencional
|
IAC10+FD
|
9.06667
|
0.15275
|
9.07 ± 0.15
|
|
Convencional
|
IAC10+H152
|
9.00000
|
0.30000
|
9 ± 0.3
|
|
Convencional
|
IAC5+FD
|
9.26667
|
0.05774
|
9.27 ± 0.06
|
|
Convencional
|
IAC5+H152
|
9.30000
|
0.17321
|
9.3 ± 0.17
|
|
Ecologico
|
BRF_PF+FD
|
9.26667
|
0.15275
|
9.27 ± 0.15
|
|
Ecologico
|
BRF_PF+H152
|
9.33333
|
0.15275
|
9.33 ± 0.15
|
|
Ecologico
|
IAC10+FD
|
9.13333
|
0.15275
|
9.13 ± 0.15
|
|
Ecologico
|
IAC10+H152
|
8.86667
|
0.05774
|
8.87 ± 0.06
|
|
Ecologico
|
IAC5+FD
|
9.40000
|
0.10000
|
9.4 ± 0.1
|
|
Ecologico
|
IAC5+H152
|
8.96667
|
0.15275
|
8.97 ± 0.15
|
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
|
Manejo
|
PortaEnxerto
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
FD
|
IAC5
|
9.40000
|
a
|
|
Ecologico
|
H152
|
BRF_PF
|
9.33333
|
a
|
|
Convencional
|
H152
|
IAC5
|
9.30000
|
ab
|
|
Convencional
|
FD
|
IAC5
|
9.26667
|
ab
|
|
Ecologico
|
FD
|
BRF_PF
|
9.26667
|
ab
|
|
Ecologico
|
FD
|
IAC10
|
9.13333
|
ab
|
|
Convencional
|
FD
|
IAC10
|
9.06667
|
ab
|
|
Convencional
|
FD
|
BRF_PF
|
9.03333
|
ab
|
|
Convencional
|
H152
|
BRF_PF
|
9.03333
|
ab
|
|
Convencional
|
H152
|
IAC10
|
9.00000
|
ab
|
|
Ecologico
|
H152
|
IAC5
|
8.96667
|
ab
|
|
Ecologico
|
H152
|
IAC10
|
8.86667
|
b
|
summary(models$models$fit_full_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0178 0.01778 0.711 0.40740
## PortaEnxerto 1 0.1111 0.11111 4.444 0.04564 *
## Copa 2 0.2956 0.14778 5.911 0.00818 **
## Manejo:PortaEnxerto 1 0.0900 0.09000 3.600 0.06987 .
## Manejo:Copa 2 0.2289 0.11444 4.578 0.02069 *
## PortaEnxerto:Copa 2 0.0956 0.04778 1.911 0.16976
## Manejo:PortaEnxerto:Copa 2 0.1067 0.05333 2.133 0.14036
## Residuals 24 0.6000 0.02500
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
9.03333
|
0.09129
|
98.95521
|
0.00000
|
|
ManejoEcologico
|
0.23333
|
0.12910
|
1.80739
|
0.08325
|
|
PortaEnxertoH152
|
0.00000
|
0.12910
|
0.00000
|
1.00000
|
|
CopaIAC10
|
0.03333
|
0.12910
|
0.25820
|
0.79845
|
|
CopaIAC5
|
0.23333
|
0.12910
|
1.80739
|
0.08325
|
|
ManejoEcologico:PortaEnxertoH152
|
0.06667
|
0.18257
|
0.36515
|
0.71820
|
|
ManejoEcologico:CopaIAC10
|
-0.16667
|
0.18257
|
-0.91287
|
0.37039
|
|
ManejoEcologico:CopaIAC5
|
-0.10000
|
0.18257
|
-0.54772
|
0.58894
|
|
PortaEnxertoH152:CopaIAC10
|
-0.06667
|
0.18257
|
-0.36515
|
0.71820
|
|
PortaEnxertoH152:CopaIAC5
|
0.03333
|
0.18257
|
0.18257
|
0.85667
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC10
|
-0.26667
|
0.25820
|
-1.03280
|
0.31200
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC5
|
-0.53333
|
0.25820
|
-2.06559
|
0.04983
|
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
BRF_PF
|
9.30000
|
a
|
|
Convencional
|
IAC5
|
9.28333
|
ab
|
|
Ecologico
|
IAC5
|
9.18333
|
ab
|
|
Convencional
|
BRF_PF
|
9.03333
|
ab
|
|
Convencional
|
IAC10
|
9.03333
|
ab
|
|
Ecologico
|
IAC10
|
9.00000
|
b
|
[[2]]
|
Manejo
|
PortaEnxerto
|
Mean
|
Group
|
|
Ecologico
|
FD
|
9.26667
|
a
|
|
Convencional
|
FD
|
9.12222
|
a
|
|
Convencional
|
H152
|
9.11111
|
a
|
|
Ecologico
|
H152
|
9.05556
|
a
|
[[3]]
|
Copa
|
PortaEnxerto
|
Mean
|
Group
|
|
IAC5
|
FD
|
9.33333
|
a
|
|
BRF_PF
|
H152
|
9.18333
|
ab
|
|
BRF_PF
|
FD
|
9.15000
|
ab
|
|
IAC5
|
H152
|
9.13333
|
ab
|
|
IAC10
|
FD
|
9.10000
|
ab
|
|
IAC10
|
H152
|
8.93333
|
b
|
summary(models$models$fit_noInt3_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0178 0.01778 0.654 0.4260
## Copa 2 0.2956 0.14778 5.437 0.0106 *
## PortaEnxerto 1 0.1111 0.11111 4.088 0.0536 .
## Manejo:Copa 2 0.2289 0.11444 4.211 0.0261 *
## Manejo:PortaEnxerto 1 0.0900 0.09000 3.311 0.0803 .
## Copa:PortaEnxerto 2 0.0956 0.04778 1.758 0.1923
## Residuals 26 0.7067 0.02718
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
8.96667
|
0.08689
|
103.19571
|
0.00000
|
|
ManejoEcologico
|
0.36667
|
0.10991
|
3.33612
|
0.00257
|
|
CopaIAC10
|
0.10000
|
0.11658
|
0.85782
|
0.39883
|
|
CopaIAC5
|
0.36667
|
0.11658
|
3.14533
|
0.00412
|
|
PortaEnxertoH152
|
0.13333
|
0.10991
|
1.21314
|
0.23598
|
|
ManejoEcologico:CopaIAC10
|
-0.30000
|
0.13461
|
-2.22867
|
0.03469
|
|
ManejoEcologico:CopaIAC5
|
-0.36667
|
0.13461
|
-2.72393
|
0.01138
|
|
ManejoEcologico:PortaEnxertoH152
|
-0.20000
|
0.10991
|
-1.81970
|
0.08034
|
|
CopaIAC10:PortaEnxertoH152
|
-0.20000
|
0.13461
|
-1.48578
|
0.14936
|
|
CopaIAC5:PortaEnxertoH152
|
-0.23333
|
0.13461
|
-1.73341
|
0.09487
|
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
BRF_PF
|
9.30000
|
a
|
|
Convencional
|
IAC5
|
9.28333
|
ab
|
|
Ecologico
|
IAC5
|
9.18333
|
ab
|
|
Convencional
|
BRF_PF
|
9.03333
|
ab
|
|
Convencional
|
IAC10
|
9.03333
|
ab
|
|
Ecologico
|
IAC10
|
9.00000
|
b
|
[[2]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
9.16111
|
a
|
|
Convencional
|
9.11667
|
a
|
[[3]]
|
Copa
|
Mean
|
Group
|
|
IAC5
|
9.23333
|
a
|
|
BRF_PF
|
9.16667
|
ab
|
|
IAC10
|
9.01667
|
b
|
[[4]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
9.19444
|
a
|
|
H152
|
9.08333
|
a
|
summary(models$models$fit_CopPort_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0178 0.01778 0.503 0.4840
## PortaEnxerto 1 0.1111 0.11111 3.142 0.0868 .
## Copa 2 0.2956 0.14778 4.179 0.0254 *
## PortaEnxerto:Copa 2 0.0956 0.04778 1.351 0.2748
## Residuals 29 1.0256 0.03536
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
9.12778
|
0.08292
|
110.07449
|
0.00000
|
|
ManejoEcologico
|
0.04444
|
0.06268
|
0.70902
|
0.48397
|
|
PortaEnxertoH152
|
0.03333
|
0.10857
|
0.30701
|
0.76103
|
|
CopaIAC10
|
-0.05000
|
0.10857
|
-0.46052
|
0.64858
|
|
CopaIAC5
|
0.18333
|
0.10857
|
1.68858
|
0.10203
|
|
PortaEnxertoH152:CopaIAC10
|
-0.20000
|
0.15354
|
-1.30255
|
0.20298
|
|
PortaEnxertoH152:CopaIAC5
|
-0.23333
|
0.15354
|
-1.51964
|
0.13943
|
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
9.16111
|
a
|
|
Convencional
|
9.11667
|
a
|
[[2]]
|
Copa
|
Mean
|
Group
|
|
IAC5
|
9.23333
|
a
|
|
BRF_PF
|
9.16667
|
ab
|
|
IAC10
|
9.01667
|
b
|
[[3]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
9.19444
|
a
|
|
H152
|
9.08333
|
a
|
summary(models$models$fit_mainEffect_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0178 0.01778 0.492 0.4885
## PortaEnxerto 1 0.1111 0.11111 3.072 0.0895 .
## Copa 2 0.2956 0.14778 4.086 0.0266 *
## Residuals 31 1.1211 0.03616
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
9.20000
|
0.07087
|
129.81073
|
0.00000
|
|
ManejoEcologico
|
0.04444
|
0.06339
|
0.70112
|
0.48846
|
|
PortaEnxertoH152
|
-0.11111
|
0.06339
|
-1.75281
|
0.08952
|
|
CopaIAC10
|
-0.15000
|
0.07764
|
-1.93207
|
0.06253
|
|
CopaIAC5
|
0.06667
|
0.07764
|
0.85870
|
0.39710
|
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))
Ratio
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Ratio")
models <- fit_qf_ccsm(dframe_qf_ccsm,'Ratio')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
Manejo
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
Convencional
|
BRF_PF+FD
|
1.56667
|
0.05774
|
1.57 ± 0.06
|
|
Convencional
|
BRF_PF+H152
|
1.96667
|
0.45092
|
1.97 ± 0.45
|
|
Convencional
|
IAC10+FD
|
1.56667
|
0.05774
|
1.57 ± 0.06
|
|
Convencional
|
IAC10+H152
|
1.93333
|
0.05774
|
1.93 ± 0.06
|
|
Convencional
|
IAC5+FD
|
1.66667
|
0.15275
|
1.67 ± 0.15
|
|
Convencional
|
IAC5+H152
|
1.76667
|
0.11547
|
1.77 ± 0.12
|
|
Ecologico
|
BRF_PF+FD
|
1.73333
|
0.05774
|
1.73 ± 0.06
|
|
Ecologico
|
BRF_PF+H152
|
1.63333
|
0.11547
|
1.63 ± 0.12
|
|
Ecologico
|
IAC10+FD
|
1.70000
|
0.26458
|
1.7 ± 0.26
|
|
Ecologico
|
IAC10+H152
|
1.56667
|
0.05774
|
1.57 ± 0.06
|
|
Ecologico
|
IAC5+FD
|
1.66667
|
0.05774
|
1.67 ± 0.06
|
|
Ecologico
|
IAC5+H152
|
1.56667
|
0.05774
|
1.57 ± 0.06
|
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
|
Manejo
|
PortaEnxerto
|
Copa
|
Mean
|
Group
|
|
Convencional
|
H152
|
BRF_PF
|
1.96667
|
a
|
|
Convencional
|
H152
|
IAC10
|
1.93333
|
a
|
|
Convencional
|
H152
|
IAC5
|
1.76667
|
a
|
|
Ecologico
|
FD
|
BRF_PF
|
1.73333
|
a
|
|
Ecologico
|
FD
|
IAC10
|
1.70000
|
a
|
|
Convencional
|
FD
|
IAC5
|
1.66667
|
a
|
|
Ecologico
|
FD
|
IAC5
|
1.66667
|
a
|
|
Ecologico
|
H152
|
BRF_PF
|
1.63333
|
a
|
|
Convencional
|
FD
|
BRF_PF
|
1.56667
|
a
|
|
Convencional
|
FD
|
IAC10
|
1.56667
|
a
|
|
Ecologico
|
H152
|
IAC10
|
1.56667
|
a
|
|
Ecologico
|
H152
|
IAC5
|
1.56667
|
a
|
summary(models$models$fit_full_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0900 0.0900 3.115 0.09028 .
## PortaEnxerto 1 0.0711 0.0711 2.462 0.12976
## Copa 2 0.0206 0.0103 0.356 0.70427
## Manejo:PortaEnxerto 1 0.3600 0.3600 12.462 0.00171 **
## Manejo:Copa 2 0.0017 0.0008 0.029 0.97160
## PortaEnxerto:Copa 2 0.0372 0.0186 0.644 0.53391
## Manejo:PortaEnxerto:Copa 2 0.0450 0.0225 0.779 0.47019
## Residuals 24 0.6933 0.0289
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.56667
|
0.09813
|
15.96511
|
0.00000
|
|
ManejoEcologico
|
0.16667
|
0.13878
|
1.20096
|
0.24149
|
|
PortaEnxertoH152
|
0.40000
|
0.13878
|
2.88231
|
0.00819
|
|
CopaIAC10
|
0.00000
|
0.13878
|
0.00000
|
1.00000
|
|
CopaIAC5
|
0.10000
|
0.13878
|
0.72058
|
0.47813
|
|
ManejoEcologico:PortaEnxertoH152
|
-0.50000
|
0.19626
|
-2.54762
|
0.01767
|
|
ManejoEcologico:CopaIAC10
|
-0.03333
|
0.19626
|
-0.16984
|
0.86656
|
|
ManejoEcologico:CopaIAC5
|
-0.16667
|
0.19626
|
-0.84921
|
0.40416
|
|
PortaEnxertoH152:CopaIAC10
|
-0.03333
|
0.19626
|
-0.16984
|
0.86656
|
|
PortaEnxertoH152:CopaIAC5
|
-0.30000
|
0.19626
|
-1.52857
|
0.13944
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC10
|
0.00000
|
0.27756
|
0.00000
|
1.00000
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC5
|
0.30000
|
0.27756
|
1.08087
|
0.29050
|
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Convencional
|
BRF_PF
|
1.76667
|
a
|
|
Convencional
|
IAC10
|
1.75000
|
a
|
|
Convencional
|
IAC5
|
1.71667
|
a
|
|
Ecologico
|
BRF_PF
|
1.68333
|
a
|
|
Ecologico
|
IAC10
|
1.63333
|
a
|
|
Ecologico
|
IAC5
|
1.61667
|
a
|
[[2]]
|
Manejo
|
PortaEnxerto
|
Mean
|
Group
|
|
Convencional
|
H152
|
1.88889
|
a
|
|
Ecologico
|
FD
|
1.70000
|
ab
|
|
Convencional
|
FD
|
1.60000
|
b
|
|
Ecologico
|
H152
|
1.58889
|
b
|
[[3]]
|
Copa
|
PortaEnxerto
|
Mean
|
Group
|
|
BRF_PF
|
H152
|
1.80000
|
a
|
|
IAC10
|
H152
|
1.75000
|
a
|
|
IAC5
|
FD
|
1.66667
|
a
|
|
IAC5
|
H152
|
1.66667
|
a
|
|
BRF_PF
|
FD
|
1.65000
|
a
|
|
IAC10
|
FD
|
1.63333
|
a
|
summary(models$models$fit_noInt3_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0900 0.0900 3.169 0.08673 .
## Copa 2 0.0206 0.0103 0.362 0.69979
## PortaEnxerto 1 0.0711 0.0711 2.504 0.12564
## Manejo:Copa 2 0.0017 0.0008 0.029 0.97111
## Manejo:PortaEnxerto 1 0.3600 0.3600 12.677 0.00145 **
## Copa:PortaEnxerto 2 0.0372 0.0186 0.655 0.52761
## Residuals 26 0.7383 0.0284
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.59167
|
0.08882
|
17.92106
|
0.00000
|
|
ManejoEcologico
|
0.11667
|
0.11234
|
1.03848
|
0.30861
|
|
CopaIAC10
|
0.00000
|
0.11916
|
0.00000
|
1.00000
|
|
CopaIAC5
|
0.02500
|
0.11916
|
0.20980
|
0.83546
|
|
PortaEnxertoH152
|
0.35000
|
0.11234
|
3.11544
|
0.00444
|
|
ManejoEcologico:CopaIAC10
|
-0.03333
|
0.13759
|
-0.24226
|
0.81048
|
|
ManejoEcologico:CopaIAC5
|
-0.01667
|
0.13759
|
-0.12113
|
0.90452
|
|
ManejoEcologico:PortaEnxertoH152
|
-0.40000
|
0.11234
|
-3.56051
|
0.00145
|
|
CopaIAC10:PortaEnxertoH152
|
-0.03333
|
0.13759
|
-0.24226
|
0.81048
|
|
CopaIAC5:PortaEnxertoH152
|
-0.15000
|
0.13759
|
-1.09018
|
0.28564
|
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Convencional
|
BRF_PF
|
1.76667
|
a
|
|
Convencional
|
IAC10
|
1.75000
|
a
|
|
Convencional
|
IAC5
|
1.71667
|
a
|
|
Ecologico
|
BRF_PF
|
1.68333
|
a
|
|
Ecologico
|
IAC10
|
1.63333
|
a
|
|
Ecologico
|
IAC5
|
1.61667
|
a
|
[[2]]
|
Manejo
|
Mean
|
Group
|
|
Convencional
|
1.74444
|
a
|
|
Ecologico
|
1.64444
|
a
|
[[3]]
|
Copa
|
Mean
|
Group
|
|
BRF_PF
|
1.72500
|
a
|
|
IAC10
|
1.69167
|
a
|
|
IAC5
|
1.66667
|
a
|
[[4]]
|
PortaEnxerto
|
Mean
|
Group
|
|
H152
|
1.73889
|
a
|
|
FD
|
1.65000
|
a
|
summary(models$models$fit_CopPort_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0900 0.09000 2.373 0.134
## PortaEnxerto 1 0.0711 0.07111 1.875 0.181
## Copa 2 0.0206 0.01028 0.271 0.765
## PortaEnxerto:Copa 2 0.0372 0.01861 0.491 0.617
## Residuals 29 1.1000 0.03793
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.70000
|
0.08588
|
19.79492
|
0.00000
|
|
ManejoEcologico
|
-0.10000
|
0.06492
|
-1.54037
|
0.13431
|
|
PortaEnxertoH152
|
0.15000
|
0.11244
|
1.33400
|
0.19259
|
|
CopaIAC10
|
-0.01667
|
0.11244
|
-0.14822
|
0.88319
|
|
CopaIAC5
|
0.01667
|
0.11244
|
0.14822
|
0.88319
|
|
PortaEnxertoH152:CopaIAC10
|
-0.03333
|
0.15902
|
-0.20962
|
0.83543
|
|
PortaEnxertoH152:CopaIAC5
|
-0.15000
|
0.15902
|
-0.94328
|
0.35333
|
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
|
Manejo
|
Mean
|
Group
|
|
Convencional
|
1.74444
|
a
|
|
Ecologico
|
1.64444
|
a
|
[[2]]
|
Copa
|
Mean
|
Group
|
|
BRF_PF
|
1.72500
|
a
|
|
IAC10
|
1.69167
|
a
|
|
IAC5
|
1.66667
|
a
|
[[3]]
|
PortaEnxerto
|
Mean
|
Group
|
|
H152
|
1.73889
|
a
|
|
FD
|
1.65000
|
a
|
summary(models$models$fit_mainEffect_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0900 0.09000 2.453 0.127
## PortaEnxerto 1 0.0711 0.07111 1.938 0.174
## Copa 2 0.0206 0.01028 0.280 0.758
## Residuals 31 1.1372 0.03668
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.73056
|
0.07138
|
24.24432
|
0.00000
|
|
ManejoEcologico
|
-0.10000
|
0.06384
|
-1.56632
|
0.12743
|
|
PortaEnxertoH152
|
0.08889
|
0.06384
|
1.39228
|
0.17375
|
|
CopaIAC10
|
-0.03333
|
0.07819
|
-0.42630
|
0.67284
|
|
CopaIAC5
|
-0.05833
|
0.07819
|
-0.74602
|
0.46128
|
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))
IT
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "IT")
models <- fit_qf_ccsm(dframe_qf_ccsm,'IT')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
Manejo
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
Convencional
|
BRF_PF+FD
|
1.84667
|
0.06807
|
1.85 ± 0.07
|
|
Convencional
|
BRF_PF+H152
|
1.80333
|
0.03215
|
1.8 ± 0.03
|
|
Convencional
|
IAC10+FD
|
1.79000
|
0.02646
|
1.79 ± 0.03
|
|
Convencional
|
IAC10+H152
|
1.73333
|
0.04726
|
1.73 ± 0.05
|
|
Convencional
|
IAC5+FD
|
1.84000
|
0.05292
|
1.84 ± 0.05
|
|
Convencional
|
IAC5+H152
|
1.92000
|
0.03606
|
1.92 ± 0.04
|
|
Ecologico
|
BRF_PF+FD
|
1.94667
|
0.07506
|
1.95 ± 0.08
|
|
Ecologico
|
BRF_PF+H152
|
1.94333
|
0.01528
|
1.94 ± 0.02
|
|
Ecologico
|
IAC10+FD
|
1.83667
|
0.03215
|
1.84 ± 0.03
|
|
Ecologico
|
IAC10+H152
|
1.64000
|
0.12124
|
1.64 ± 0.12
|
|
Ecologico
|
IAC5+FD
|
1.95333
|
0.04163
|
1.95 ± 0.04
|
|
Ecologico
|
IAC5+H152
|
1.79667
|
0.09018
|
1.8 ± 0.09
|
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
|
Manejo
|
PortaEnxerto
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
FD
|
IAC5
|
1.95333
|
a
|
|
Ecologico
|
FD
|
BRF_PF
|
1.94667
|
a
|
|
Ecologico
|
H152
|
BRF_PF
|
1.94333
|
a
|
|
Convencional
|
H152
|
IAC5
|
1.92000
|
a
|
|
Convencional
|
FD
|
BRF_PF
|
1.84667
|
ab
|
|
Convencional
|
FD
|
IAC5
|
1.84000
|
ab
|
|
Ecologico
|
FD
|
IAC10
|
1.83667
|
ab
|
|
Convencional
|
H152
|
BRF_PF
|
1.80333
|
abc
|
|
Ecologico
|
H152
|
IAC5
|
1.79667
|
abc
|
|
Convencional
|
FD
|
IAC10
|
1.79000
|
abc
|
|
Convencional
|
H152
|
IAC10
|
1.73333
|
bc
|
|
Ecologico
|
H152
|
IAC10
|
1.64000
|
c
|
summary(models$models$fit_full_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.00840 0.00840 2.283 0.14385
## PortaEnxerto 1 0.03547 0.03547 9.637 0.00484 **
## Copa 2 0.13815 0.06908 18.768 1.24e-05 ***
## Manejo:PortaEnxerto 1 0.02834 0.02834 7.699 0.01053 *
## Manejo:Copa 2 0.03651 0.01825 4.959 0.01575 *
## PortaEnxerto:Copa 2 0.01871 0.00935 2.541 0.09977 .
## Manejo:PortaEnxerto:Copa 2 0.02957 0.01479 4.017 0.03127 *
## Residuals 24 0.08833 0.00368
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.84667
|
0.03503
|
52.72207
|
0.00000
|
|
ManejoEcologico
|
0.10000
|
0.04953
|
2.01878
|
0.05482
|
|
PortaEnxertoH152
|
-0.04333
|
0.04953
|
-0.87480
|
0.39035
|
|
CopaIAC10
|
-0.05667
|
0.04953
|
-1.14398
|
0.26392
|
|
CopaIAC5
|
-0.00667
|
0.04953
|
-0.13459
|
0.89406
|
|
ManejoEcologico:PortaEnxertoH152
|
0.04000
|
0.07005
|
0.57100
|
0.57331
|
|
ManejoEcologico:CopaIAC10
|
-0.05333
|
0.07005
|
-0.76133
|
0.45388
|
|
ManejoEcologico:CopaIAC5
|
0.01333
|
0.07005
|
0.19033
|
0.85065
|
|
PortaEnxertoH152:CopaIAC10
|
-0.01333
|
0.07005
|
-0.19033
|
0.85065
|
|
PortaEnxertoH152:CopaIAC5
|
0.12333
|
0.07005
|
1.76057
|
0.09105
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC10
|
-0.18000
|
0.09907
|
-1.81690
|
0.08174
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC5
|
-0.27667
|
0.09907
|
-2.79265
|
0.01010
|
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
BRF_PF
|
1.94500
|
a
|
|
Convencional
|
IAC5
|
1.88000
|
ab
|
|
Ecologico
|
IAC5
|
1.87500
|
ab
|
|
Convencional
|
BRF_PF
|
1.82500
|
bc
|
|
Convencional
|
IAC10
|
1.76167
|
bc
|
|
Ecologico
|
IAC10
|
1.73833
|
c
|
[[2]]
|
Manejo
|
PortaEnxerto
|
Mean
|
Group
|
|
Ecologico
|
FD
|
1.91222
|
a
|
|
Convencional
|
FD
|
1.82556
|
ab
|
|
Convencional
|
H152
|
1.81889
|
b
|
|
Ecologico
|
H152
|
1.79333
|
b
|
[[3]]
|
Copa
|
PortaEnxerto
|
Mean
|
Group
|
|
BRF_PF
|
FD
|
1.89667
|
a
|
|
IAC5
|
FD
|
1.89667
|
a
|
|
BRF_PF
|
H152
|
1.87333
|
a
|
|
IAC5
|
H152
|
1.85833
|
a
|
|
IAC10
|
FD
|
1.81333
|
a
|
|
IAC10
|
H152
|
1.68667
|
b
|
summary(models$models$fit_noInt3_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.00840 0.00840 1.853 0.18512
## Copa 2 0.13815 0.06908 15.232 4.18e-05 ***
## PortaEnxerto 1 0.03547 0.03547 7.822 0.00958 **
## Manejo:Copa 2 0.03651 0.01825 4.025 0.03000 *
## Manejo:PortaEnxerto 1 0.02834 0.02834 6.249 0.01907 *
## Copa:PortaEnxerto 2 0.01871 0.00935 2.062 0.14744
## Residuals 26 0.11791 0.00453
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.80861
|
0.03549
|
50.95842
|
0.00000
|
|
ManejoEcologico
|
0.17611
|
0.04489
|
3.92281
|
0.00057
|
|
CopaIAC10
|
-0.01167
|
0.04762
|
-0.24501
|
0.80837
|
|
CopaIAC5
|
0.06250
|
0.04762
|
1.31255
|
0.20080
|
|
PortaEnxertoH152
|
0.03278
|
0.04489
|
0.73011
|
0.47185
|
|
ManejoEcologico:CopaIAC10
|
-0.14333
|
0.05498
|
-2.60683
|
0.01493
|
|
ManejoEcologico:CopaIAC5
|
-0.12500
|
0.05498
|
-2.27340
|
0.03150
|
|
ManejoEcologico:PortaEnxertoH152
|
-0.11222
|
0.04489
|
-2.49971
|
0.01907
|
|
CopaIAC10:PortaEnxertoH152
|
-0.10333
|
0.05498
|
-1.87934
|
0.07145
|
|
CopaIAC5:PortaEnxertoH152
|
-0.01500
|
0.05498
|
-0.27281
|
0.78716
|
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Ecologico
|
BRF_PF
|
1.94500
|
a
|
|
Convencional
|
IAC5
|
1.88000
|
a
|
|
Ecologico
|
IAC5
|
1.87500
|
a
|
|
Convencional
|
BRF_PF
|
1.82500
|
a
|
|
Convencional
|
IAC10
|
1.76167
|
ab
|
|
Ecologico
|
IAC10
|
1.73833
|
b
|
[[2]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
1.85278
|
a
|
|
Convencional
|
1.82222
|
a
|
[[3]]
|
Copa
|
Mean
|
Group
|
|
BRF_PF
|
1.8850
|
a
|
|
IAC5
|
1.8775
|
a
|
|
IAC10
|
1.7500
|
b
|
[[4]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
1.86889
|
a
|
|
H152
|
1.80611
|
b
|
summary(models$models$fit_CopPort_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.00840 0.00840 1.333 0.257623
## PortaEnxerto 1 0.03547 0.03547 5.629 0.024519 *
## Copa 2 0.13815 0.06908 10.961 0.000285 ***
## PortaEnxerto:Copa 2 0.01871 0.00935 1.484 0.243402
## Residuals 29 0.18275 0.00630
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.88139
|
0.03500
|
53.74702
|
0.00000
|
|
ManejoEcologico
|
0.03056
|
0.02646
|
1.15474
|
0.25762
|
|
PortaEnxertoH152
|
-0.02333
|
0.04583
|
-0.50911
|
0.61453
|
|
CopaIAC10
|
-0.08333
|
0.04583
|
-1.81825
|
0.07937
|
|
CopaIAC5
|
0.00000
|
0.04583
|
0.00000
|
1.00000
|
|
PortaEnxertoH152:CopaIAC10
|
-0.10333
|
0.06482
|
-1.59426
|
0.12172
|
|
PortaEnxertoH152:CopaIAC5
|
-0.01500
|
0.06482
|
-0.23143
|
0.81861
|
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
|
Manejo
|
Mean
|
Group
|
|
Ecologico
|
1.85278
|
a
|
|
Convencional
|
1.82222
|
a
|
[[2]]
|
Copa
|
Mean
|
Group
|
|
BRF_PF
|
1.8850
|
a
|
|
IAC5
|
1.8775
|
a
|
|
IAC10
|
1.7500
|
b
|
[[3]]
|
PortaEnxerto
|
Mean
|
Group
|
|
FD
|
1.86889
|
a
|
|
H152
|
1.80611
|
b
|
summary(models$models$fit_mainEffect_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.00840 0.00840 1.293 0.264202
## PortaEnxerto 1 0.03547 0.03547 5.458 0.026119 *
## Copa 2 0.13815 0.06908 10.629 0.000305 ***
## Residuals 31 0.20145 0.00650
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
1.90111
|
0.03004
|
63.28022
|
0.00000
|
|
ManejoEcologico
|
0.03056
|
0.02687
|
1.13712
|
0.26420
|
|
PortaEnxertoH152
|
-0.06278
|
0.02687
|
-2.33626
|
0.02612
|
|
CopaIAC10
|
-0.13500
|
0.03291
|
-4.10208
|
0.00027
|
|
CopaIAC5
|
-0.00750
|
0.03291
|
-0.22789
|
0.82123
|
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))
VitC
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "VitC")
models <- fit_qf_ccsm(dframe_qf_ccsm,'VitC')
Graficos
Barplot 1
results$barp

Barplot 2
results$barp_2

Barplot + SD 1
results$barp_error

Barplot + SD 2
results$barp_error_2

Mean Error 1
results$mean_desv_plot

Mean Error 2
results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
|
Manejo
|
Tratamento
|
Mean
|
SD
|
Mean ± SD
|
|
Convencional
|
BRF_PF+FD
|
34.15667
|
2.19805
|
34.16 ± 2.2
|
|
Convencional
|
BRF_PF+H152
|
37.52333
|
1.58418
|
37.52 ± 1.58
|
|
Convencional
|
IAC10+FD
|
34.25333
|
1.33121
|
34.25 ± 1.33
|
|
Convencional
|
IAC10+H152
|
33.71000
|
2.57874
|
33.71 ± 2.58
|
|
Convencional
|
IAC5+FD
|
32.83667
|
2.87114
|
32.84 ± 2.87
|
|
Convencional
|
IAC5+H152
|
33.98333
|
1.23678
|
33.98 ± 1.24
|
|
Ecologico
|
BRF_PF+FD
|
31.96000
|
1.10585
|
31.96 ± 1.11
|
|
Ecologico
|
BRF_PF+H152
|
33.32333
|
3.61777
|
33.32 ± 3.62
|
|
Ecologico
|
IAC10+FD
|
34.23000
|
2.17690
|
34.23 ± 2.18
|
|
Ecologico
|
IAC10+H152
|
36.66667
|
0.76846
|
36.67 ± 0.77
|
|
Ecologico
|
IAC5+FD
|
31.53667
|
1.24340
|
31.54 ± 1.24
|
|
Ecologico
|
IAC5+H152
|
33.23333
|
1.44901
|
33.23 ± 1.45
|
Modelos
Full
models$tukey$fit_full_aov %>% kable_cust()
|
Manejo
|
PortaEnxerto
|
Copa
|
Mean
|
Group
|
|
Convencional
|
H152
|
BRF_PF
|
37.52333
|
a
|
|
Ecologico
|
H152
|
IAC10
|
36.66667
|
ab
|
|
Convencional
|
FD
|
IAC10
|
34.25333
|
ab
|
|
Ecologico
|
FD
|
IAC10
|
34.23000
|
ab
|
|
Convencional
|
FD
|
BRF_PF
|
34.15667
|
ab
|
|
Convencional
|
H152
|
IAC5
|
33.98333
|
ab
|
|
Convencional
|
H152
|
IAC10
|
33.71000
|
ab
|
|
Ecologico
|
H152
|
BRF_PF
|
33.32333
|
ab
|
|
Ecologico
|
H152
|
IAC5
|
33.23333
|
ab
|
|
Convencional
|
FD
|
IAC5
|
32.83667
|
ab
|
|
Ecologico
|
FD
|
BRF_PF
|
31.96000
|
ab
|
|
Ecologico
|
FD
|
IAC5
|
31.53667
|
b
|
summary(models$models$fit_full_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 7.60 7.599 1.868 0.1844
## PortaEnxerto 1 22.40 22.404 5.507 0.0275 *
## Copa 2 21.33 10.665 2.621 0.0934 .
## Manejo:PortaEnxerto 1 0.58 0.583 0.143 0.7084
## Manejo:Copa 2 32.69 16.347 4.018 0.0313 *
## PortaEnxerto:Copa 2 3.13 1.564 0.384 0.6850
## Manejo:PortaEnxerto:Copa 2 9.31 4.657 1.145 0.3351
## Residuals 24 97.64 4.069
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_full)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
34.15667
|
1.16455
|
29.33034
|
0.00000
|
|
ManejoEcologico
|
-2.19667
|
1.64692
|
-1.33380
|
0.19479
|
|
PortaEnxertoH152
|
3.36667
|
1.64692
|
2.04422
|
0.05205
|
|
CopaIAC10
|
0.09667
|
1.64692
|
0.05870
|
0.95368
|
|
CopaIAC5
|
-1.32000
|
1.64692
|
-0.80149
|
0.43071
|
|
ManejoEcologico:PortaEnxertoH152
|
-2.00333
|
2.32910
|
-0.86013
|
0.39823
|
|
ManejoEcologico:CopaIAC10
|
2.17333
|
2.32910
|
0.93312
|
0.36005
|
|
ManejoEcologico:CopaIAC5
|
0.89667
|
2.32910
|
0.38498
|
0.70364
|
|
PortaEnxertoH152:CopaIAC10
|
-3.91000
|
2.32910
|
-1.67876
|
0.10617
|
|
PortaEnxertoH152:CopaIAC5
|
-2.22000
|
2.32910
|
-0.95316
|
0.35001
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC10
|
4.98333
|
3.29385
|
1.51292
|
0.14336
|
|
ManejoEcologico:PortaEnxertoH152:CopaIAC5
|
2.55333
|
3.29385
|
0.77518
|
0.44580
|
par(mfrow=c(2,2))
plot(models$models$fit_full_aov)

par(mfrow=c(1,1))
2x2 Interaction
lapply(models$tukey$fit_noInt3_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Convencional
|
BRF_PF
|
35.84000
|
a
|
|
Ecologico
|
IAC10
|
35.44833
|
a
|
|
Convencional
|
IAC10
|
33.98167
|
a
|
|
Convencional
|
IAC5
|
33.41000
|
a
|
|
Ecologico
|
BRF_PF
|
32.64167
|
a
|
|
Ecologico
|
IAC5
|
32.38500
|
a
|
[[2]]
|
Manejo
|
PortaEnxerto
|
Mean
|
Group
|
|
Convencional
|
H152
|
35.07222
|
a
|
|
Ecologico
|
H152
|
34.40778
|
a
|
|
Convencional
|
FD
|
33.74889
|
a
|
|
Ecologico
|
FD
|
32.57556
|
a
|
[[3]]
|
Copa
|
PortaEnxerto
|
Mean
|
Group
|
|
BRF_PF
|
H152
|
35.42333
|
a
|
|
IAC10
|
H152
|
35.18833
|
a
|
|
IAC10
|
FD
|
34.24167
|
a
|
|
IAC5
|
H152
|
33.60833
|
a
|
|
BRF_PF
|
FD
|
33.05833
|
a
|
|
IAC5
|
FD
|
32.18667
|
a
|
summary(models$models$fit_noInt3_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 7.60 7.599 1.847 0.1858
## Copa 2 21.33 10.665 2.593 0.0940 .
## PortaEnxerto 1 22.40 22.404 5.446 0.0276 *
## Manejo:Copa 2 32.69 16.347 3.974 0.0312 *
## Manejo:PortaEnxerto 1 0.58 0.583 0.142 0.7097
## Copa:PortaEnxerto 2 3.13 1.564 0.380 0.6875
## Residuals 26 106.96 4.114
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_noInt3)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
34.78472
|
1.06898
|
32.53997
|
0.00000
|
|
ManejoEcologico
|
-3.45278
|
1.35217
|
-2.55351
|
0.01687
|
|
CopaIAC10
|
-1.14917
|
1.43419
|
-0.80126
|
0.43024
|
|
CopaIAC5
|
-1.95833
|
1.43419
|
-1.36546
|
0.18381
|
|
PortaEnxertoH152
|
2.11056
|
1.35217
|
1.56087
|
0.13065
|
|
ManejoEcologico:CopaIAC10
|
4.66500
|
1.65606
|
2.81692
|
0.00914
|
|
ManejoEcologico:CopaIAC5
|
2.17333
|
1.65606
|
1.31235
|
0.20087
|
|
ManejoEcologico:PortaEnxertoH152
|
0.50889
|
1.35217
|
0.37635
|
0.70971
|
|
CopaIAC10:PortaEnxertoH152
|
-1.41833
|
1.65606
|
-0.85645
|
0.39957
|
|
CopaIAC5:PortaEnxertoH152
|
-0.94333
|
1.65606
|
-0.56962
|
0.57382
|
par(mfrow=c(2,2))
plot(models$models$fit_noInt3_aov)

par(mfrow=c(1,1))
Copa*PortaEnxerto Interaction + Manejo
lapply(models$tukey$fit_CopPort_aov,kable_cust)
[[1]]
|
Manejo
|
Copa
|
Mean
|
Group
|
|
Convencional
|
BRF_PF
|
35.84000
|
a
|
|
Ecologico
|
IAC10
|
35.44833
|
a
|
|
Convencional
|
IAC10
|
33.98167
|
a
|
|
Convencional
|
IAC5
|
33.41000
|
a
|
|
Ecologico
|
BRF_PF
|
32.64167
|
a
|
|
Ecologico
|
IAC5
|
32.38500
|
a
|
[[2]]
|
Manejo
|
Mean
|
Group
|
|
Convencional
|
34.41056
|
a
|
|
Ecologico
|
33.49167
|
a
|
[[3]]
|
Copa
|
Mean
|
Group
|
|
IAC10
|
34.71500
|
a
|
|
BRF_PF
|
34.24083
|
a
|
|
IAC5
|
32.89750
|
a
|
[[4]]
|
PortaEnxerto
|
Mean
|
Group
|
|
H152
|
34.74000
|
a
|
|
FD
|
33.16222
|
b
|
summary(models$models$fit_CopPort_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 7.60 7.599 1.571 0.2200
## PortaEnxerto 1 22.40 22.404 4.633 0.0398 *
## Copa 2 21.33 10.665 2.206 0.1283
## PortaEnxerto:Copa 2 3.13 1.564 0.323 0.7263
## Residuals 29 140.24 4.836
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_CopPort)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
33.51778
|
0.96968
|
34.56581
|
0.00000
|
|
ManejoEcologico
|
-0.91889
|
0.73301
|
-1.25358
|
0.22001
|
|
PortaEnxertoH152
|
2.36500
|
1.26961
|
1.86278
|
0.07265
|
|
CopaIAC10
|
1.18333
|
1.26961
|
0.93205
|
0.35901
|
|
CopaIAC5
|
-0.87167
|
1.26961
|
-0.68656
|
0.49781
|
|
PortaEnxertoH152:CopaIAC10
|
-1.41833
|
1.79550
|
-0.78994
|
0.43598
|
|
PortaEnxertoH152:CopaIAC5
|
-0.94333
|
1.79550
|
-0.52539
|
0.60331
|
par(mfrow=c(2,2))
plot(models$models$fit_CopPort_aov)

par(mfrow=c(1,1))
Main Effects
lapply(models$tukey$fit_mainEffect_aov,kable_cust)
[[1]]
|
Manejo
|
Mean
|
Group
|
|
Convencional
|
34.41056
|
a
|
|
Ecologico
|
33.49167
|
a
|
[[2]]
|
Copa
|
Mean
|
Group
|
|
IAC10
|
34.71500
|
a
|
|
BRF_PF
|
34.24083
|
a
|
|
IAC5
|
32.89750
|
a
|
[[3]]
|
PortaEnxerto
|
Mean
|
Group
|
|
H152
|
34.74000
|
a
|
|
FD
|
33.16222
|
b
|
summary(models$models$fit_mainEffect_aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 7.60 7.599 1.643 0.2094
## PortaEnxerto 1 22.40 22.404 4.845 0.0353 *
## Copa 2 21.33 10.665 2.306 0.1165
## Residuals 31 143.36 4.625
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(models$models$fit_mainEffect)$coeff %>% kable_cust()
|
|
Estimate
|
Std. Error
|
t value
|
Pr(>|t|)
|
|
(Intercept)
|
33.91139
|
0.80144
|
42.31301
|
0.00000
|
|
ManejoEcologico
|
-0.91889
|
0.71683
|
-1.28188
|
0.20939
|
|
PortaEnxertoH152
|
1.57778
|
0.71683
|
2.20105
|
0.03531
|
|
CopaIAC10
|
0.47417
|
0.87793
|
0.54009
|
0.59299
|
|
CopaIAC5
|
-1.34333
|
0.87793
|
-1.53011
|
0.13613
|
par(mfrow=c(2,2))
plot(models$models$fit_mainEffect_aov)

par(mfrow=c(1,1))